[ 251 ] 


VII. The Stability of the Pear-Shaped Figure of Equilibrium of a Rotating 

Mass of Liquid . 

By G. H. Darwin, F.R.S ., Plumian Professor and Fellow of Trinity College , 

in the University of Cambridge. 

Received June 19,—Read June 19, 1902, 


Introduction 


Table of Contents, 


Pago 

252 


Part I.—Analytical Investigation, 


§ 1. Method of procedure. . ..•.. 253 

§ 2. The lost energy of the system. ., . . . . . .. . . 254 

§ 3, Expression for the element of volume ... . . , ... 256 

§ 4. Determination of h; definition of symbols for integrals ....... 260 

§ 5. The energies \JJ and JR . 262 

§ 6. Surface density of concentration G; energy GR . . . . . . . . . 264 

§ 7. The energy lGO; result for \JJ - JR + GR - \GG ... . . . . . 267 

§ 8. The term -^Mdho 2 ....268 

§ 9. Double layers.... 268 

§ 10. Determination of e and X.274 

§11. The energy §7rp 2 Je 3 (1 -- Xe) d<r ...275 

* 7 77" 

§12. The energy \p e 2 — dm. ... . 276 

§13. Terms in the energy depending on the moment of inertia. , . . . . . 279 

§ 14. The lost energy of the system; solution of the problem . , . . , . . 282 


Part II. —Numerical Calculation. 


§ 15. Determination of certain integrals . . . ... 284 

§ 16. The integrals <r 2 , <r 4 , ... ... . 290 

§ 17. The integrals w/, pf, .. 291 

§ 18. The integrals ...294 

§19. Synthesis of numerical results; stability of the pear . . ... . . . 296 

§ 20. Second approximation to the form of the pear.. 301 

Part III. 

Summary. ...306 

(327.) 2 K 2 9-1.03 


The Royal Society is collaborating with JSTOR to digitize, preserve, and extend access to 
Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character. 

www.jstor.org 





























252 PROFESSOR G, H, DARWIN ON THE STABILITY OF THE PEAR-SHAPED 


Introduction. 

By aid of the methods of a paper on “Ellipsoidal Harmonic Analysis” (‘Phil. Trans., 
A, vol. 197, pp. 461-557), I here resume the subject of a previous paper (‘ Phil. Trans./ 
A, vol. 198, pp. 301-831). These papers will lie referred to hereafter by the abridged 
titles of “ Harmonics ” and “ The Pear-shaped Figure.” 

At the end of the latter of these it was stated that the stability of the figure could 
not be proved definitely without approximation of a higher order of accuracy. A fter 
some correspondence with M. Poincare during the course of my work, I made an 
attempt to carry out this further approximation, but found that the expression 
for a certain portion of the energy entirely foiled me. Meanwhile he had turned his 
attention to the subject, and he has shown (‘Phil. Trans.,’A, vol. 198, pp. 333-373) 
by a method of the greatest ingenuity and skill how the problem may be solved. He 
has not, however, pursued the arduous task of converting his analytical results into 
numbers, so that he left the question as to the stability of the pear still unanswered. 

M. Poincare was so kind as to allow me to detain his manuscript on its way to the 
'Royal Society tor two or three days, and I devoted that time almost entirely to 
understanding the method of his attack on the key of the position—namely, the 
method of double layers, expounded in my own language in § 9 below. Being thus 
furnished with the means, I was able to resume my attempt under favourable 
conditions, and this paper is the result. 

The substance of the analysis of this paper is, of course, essentially the same as his, 
but the arrangement and notation are so different that the two present but little 
superficial resemblance. This difference arises partly from the fact that I desired to 
use my own notation for the ellipsoidal harmonics, and partly because during the time 
that I was working at the analysis his paper was still unprinted and therefore 
inaccessible to me. But it is, perhaps, well that the two investigations of so 
complicated a subject should be nearly independent of one another. 

It is rather unfortunate that I did not feel myself sufficiently expert in the use of 
the methods of Weierstrass and Schwarz to evaluate the elliptic integrals after the 
methods suggested by M. Poincare, but every exertion has been taken to insure 
correctness in the arithmetical results, on. which the proof of stability depends. My 
choice of antiquated methods of computation leaves the way open for some one else to 
verify the conclusions by wholly independent and more elegant calculations. It is 
highly desirable that such a verification should be made. 

As the body of this paper will hardly be studied by any one unless they should he 
actually working at the subject, I give a summary at the end. Even the mathe¬ 
matician who desires to study the subject in detail may find it advantageous to read 
the summary before looking at the analytical investigation. 
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PART I. 

Analytical Investigation. 

§ I. Method of Procedure. 

The pear-shaped figure is a deformation ot the critical Jacobian ellipsoid, and to 
the first order of small quantities it is expressed by the third zonal harmonic with 
respect to the longest axis of the ellipsoid. In the higher approximation a number 
of other harmonic terms will arise, and the coefficients of these new terms will be ot 
the second order of small quantities. The mass of an harmonic inequality vanishes 
only to the first order, and it can no longer be assumed that the centre of inertia of 
the pear coincides with the centre of the ellipsoid. 

In order to define the pear, I describe an ellipsoid similar to and concentric with 
the original critical Jacobian; this new ellipsoid is taken to be sufficiently large to 
enclose the whole of the pear. It is clearly itself a critical Jacobian, and I adopt it 
as the ellipsoid of reference, and call it J. I call the region between J and the 
pear R . The pear may then be defined by density -J- p throughout J, and density — p 
throughout R. 

If h is the parameter which defines ./, its axes are expressed in the notation of 
“ Harmonics” by Jcv 0 , k (yf — 1) J , k [vf — - - —j ; or in the notation of the “ Pear- 

shaped Figure ” by /c/sin /3, Jc cos /3/sin /3, h cos y/sin /3, where sin /3 = k sin y. 

Now let Si a denote any surface harmonic, so that S? is the same thing as 
DPd (/ x ) or P i s (f-)] X [Cj (<f>) or Cj (</>)]. The third zonal harmonic deformation will 
then be eS s or eflj (p) C 3 (</>), where e is of the first order of small quantities. On 
account of the symmetry of the figure, the new terms cannot involve the sine 
functions % or S, and moreover, the rank s must necessarily be even. 

Suppose that the new terms are expressed by %ffSf for all values of i from 1 to 
infinity, and with s equal to 0, 2, 4 . . . i or i — 1. Then all the f s are of order e 3 , 
excepting jf> which is zero. 

We have seen in “ Harmonics,” §11, that if p 0 denotes the perpendicular from the 
centre of the ellipsoid v 0 on to the tangent plane at /x, </>, the equation to a harmonic 
deformation of the ellipsoid is 


(f' — I'o 2 ) = 2&SJ. 

AO 

Since this equation may be written in the form 


O _ 1 t & W - t) ® 

r 1 - pj 


+ 70 o — 1 + 2 eS; s 
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it is clear that if 2eSf is a constant, say c, the surface defined is an ellipsoid similar 
to the surface of reference, with semi-axes augmented in the proportion of (1 + c) i: to 
unity. 

I now replace the variable v by a new one, namely, 


T 


& 

w 




2 \ 

v q) 


( 1 ). 


The negative sign is taken because the points to he specified will lie inside J. 

Then r = c, a constant, defines an interior ellipsoid similar to and concentric with J. 
The equation to the pear may now be written 

t = c — eS ?) — '£fi s S; s . 


Th only condition which has been imposed on c is that it shall be great enough to 
make r always positive. 

In order to solve our problem it is necessary to determine the energy lost in the 
process of concentration from a condition of infinite dispersion into the final con¬ 
figuration. This involves the use of the formula for the gravity of J, inclusive of 
rotation. It is well known that this formula is simple for the inside of J and more 
complicated for the outside. Since the whole region R lies inside J there is no 
necessity in the present case to use the more complicated formula. 

The final expression for the lost energy cannot involve the size of J, the exterior 
ellipsoid of reference, and therefore the arbitrary constant c must ultimately disappear. 
It is therefore legitimate to make c zero from the beginning. 

It is clear that we might with equal justice have discussed the problem by means 
of an ellipsoid which should lie entirely inside the pear, the region between the pear 
and the ellipsoid would then have been filled with positive density, and the formula 
for external gravity would have been needed. The same argument as before would 
then have justified our putting the constant c equal to zero. 

We thus arrive at the same conclusion as does M. Poincare, namely, that it is 
immaterial whether the formula for external or internal gravity be used. 

I now revert to my first hypothesis of the enveloping ellipsoid, but put c equal 
to zero from the first. In order, however, to afford clearness to our conceptions, I 
shall continue to discuss the problem as though c were not zero and as though J 
enclosed the whole pear. With this explanation, we may write the equation to the 
pear in the form 

oo 

t = — eS 3 — Zf'Sf ..(2). 


§ 2. The Lost Energy of the System. 

If the negative density in R is transported along tubes formed by a family ot 
orthogonal curves and deposited as surface density on J, we may refer to such a 
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condensation as — C. I do not suppose the condensation actually effected, but 
imagine the surface of J to be coated with equal and opposite condensations + G 
and — C. 

The system of masses forming the pear may then be considered as being as 
follows:— 

Density + p throughout J, say + J. 

Negative condensation on J, say — C.. 

Positive condensation + C on J and negative volume density — p throughout R. 
This last forms a double system of zero mass, say I), and D — C — R. 

Let Vj, V r be the potentials of + J and + R, and the potential of the pear. 

An element of volume being written dv, let dv, do, dv denote integrations 

j Jr Jj-r 

throughout J, R and the pear respectively. 

Let d be the distance along the z axis from the centre of the ellipsoid as origin to 
the centre of inertia of the pear; let a> be the angular velocity of the critical 
Jacobian about the axis x, so that co 3 / 27 rp = T4200; and let <w 3 + §w 3 be the square 
of the angular velocity of the pear. Lastly, let M be the mass of the pear. 

Then the lost energy E is given by 


E=d Vj_ r pdv + | (or + So, 3 ) f [ ?/ 3 + (z - df\ pdv. 
Jj->' Jj-r 

Now ( zpdv = Md, so that f (— 2 zd + #) pdv — — MR. 

Jj- r Jj~> 


Again, since 


we have 


Vj-r = Vj - V n f = f. - f , f. V rp dv = f Vjpdv, 

J j—r Jj Jr J j Jr 

\ Vj_ r pdv = | Vjpdv — Vjpdv + V r pdv. 

J j—r Jj Jr Jr 


Also 


Jr (<u 3 + 8 oj 3 ) | [?/ 3 -}- (z — o?) 3 ] pdv — |&) 3 j ( y 3 + « 3 ) pdv — ^w 3 j" (t/ 3 + z z ) pdv 

+ t|S w 3 f (y 3 + * 3 ) pdv - i (®8 + BR) Md\ 
Jj-r 

Hence 

E = 4 ( . [ Vj + «, 8 (f + **)] pdv - [ [ Vj + (/ +#)] pdv + i £ Vrpdv 

+ 48<» 8 f (f + z s ) pdv - I (of + hvf) Md 8 

Jj-r 


As the several terms will be considered separately, it will be convenient to have an 
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abridged notation to specify them. I may denote the lost energy of J, inclusive of 
rotation, hy \JJ ; the mutual lost energy of J and of the region R, considered as 
filled with positive density, by JR ; the lost energy of the region R by },RR. 

The moment of inertia of the pear is A, and it is equal to Aj — A r , the moment of 
inertia of J less that of R. 

Then 

E = \JJ - JR + %RR + J (Aj - A,) hcJ - i (oJ + So> 2 ) Md\ 

where 

i«n = H [Pi-Wfe/ 8 + **)]>>*, 

J j 

JR — f W + W (*,* + #)] pdv, 

J r 

Aj — [ (if + z 2 ) pdv, A r — fy 3 + z~) pdv, 

and 

fRR = f V,pdv. 

If \DD denotes the lost energy of the double system described above, we clearly 
have 

\RR = i(C- R)(C- R) + OR - \GG = fDD + OR - f CC. 

We require to evaluate E to the fourth order; now d is at least of the second 
order and d 3 of the fourth order; hence d 3 . 8co' 2 is at least of the fifth order and 
negligible. 

Hence, finally, to the required degree of approximation 

E = \JJ - JR + GR - ±CC + \I)D + \ (A, - A r ) - \Md?<J . . (3). 

It will appear below that d is not even of the second order, so that the last term 
will, in fact, entirely disappear, although we cannot see at the present stage that this 
will he so. 


3. Expression for the Element of Volume. 


The parameter 0 of “ Harmonics ” is connected with k of the “ Pear-sha 
hy the equations 


■'lire 


1-0 a 1 -k 2 J 2 20 J 2 2/3 Jjt 

1 + 0 K ’ P - 1 + *2 “ 1. + ** ‘ T-0 ~~ 1 + 0 “ * ' 

There will, I think, be no confusion if I also use 0 in a second sense, defining it by 
the equations 


sin 0 = k sin y, 


cos 3 0=1 — k~ sin® y. 
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It has already been remarked above that the squares of the semi-axes of J are 


Z; 2 

*v = -. 

u sin 2 p 


k 2 cos 2 /3 
sin 3 /3 9 




1 4- 6\ _ Id cos 3 7 

1-/3/ sin 3 /3 


The mass of J is then | npk 


3 cos /3 cos 7 


sin 3 /3 


I now take the mass M of the pear to be 

i>r a 7 , cos /3 cos 7 
M = inpkJ —: , - \ 
d ' 0 sin 8 /3 

Thus k 0 is a constant which specifies the volume of liquid in the pear, and the mass 
of JisM (k/k 0 f. 

It will be convenient to introduce certain new symbols, namely, 

A 3 = 1 — k 3 sin 3 6, r 3 = 1 — k / 3 cos 3 <j>, 

A a 3 = 1 — k 3 sin 3 y sin 3 8, W = cos 3 y + k h sin 3 y cos 3 (f>, 

where sin 8 is the /x of “ Harmonics.” 

The roots of the fundamental cubic were n 3 , p 3 , and and in the new 


1 — K n ~ cos 3 d> r 2 

notation they are r 3 , sin 3 8, -~— s ~ f-— or 2 . 

fC n ' K 


sin 3 /3 5 


Since zy 3 = . 0 , we now have 


v o “ K 


.A’ 

sin 3 /S ’ 


1-/3 cos 20_I\ 3 

1 — /3 sin 2 /3 ’ 


The expression for p 0 , the perpendicular from the centre on to the tangent plane 
at 6, (j), is given in (49) of “ Harmonics,” namely, 


Pd 

Id 


vw-u (*-££) 

1 — j3 cos 2<£ 


W - fQ r 0 2 - 


1-/8 


cos 2 j3 cos 2 7 1 

" : sin 2 /9 # * * 


(4). 


Also by (50) of “ Harmonics” the element of surface da of the ellipsoid is given by 

1 — (3 cos 2(/> 


?JL da _ /.: 


d<f> 


k\ (n 0 3 — 1)* v, 


1 + £\* 


1-/8 


- p- 


1-/3/ /1 — /3 cos 2 (fr A / I + 8 _ o 

1-/3 / U-/3. * 


Passing to the new notation this may be written 

* \ 8 1 - * 2 s in 2 6 - k' 2 cos 2 j> _ 3Jf / k y A^ U, 2 / I 


d6 d(f> 4:irp \ Jc 0 
YOL. CC.—A. 


AU 


4:7rp \ /% / * AU sin 2 7 Viy Aj 2 / * 


2 L 
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The new independent variable t is to replace v ; it was defined in (1) by 

It? 


2p ( 


w —*?)> 


and in accordance with (2) the equation to the surface of the pear is 

r = — eS s — t ffSf. 

1 

From (4) 

a 9 2p 3 1 / 2 t cos 2 8 cos 3 y\ 

For brevity I now write 


V1V 


2 t cos 2 8 cos 2 7 . i . 

--:— . SO that 


A^iy 


sin 2 /3 




sin 2 /3 


A, 2 


1 


sin 2 /3 \ Af 


1 — /3 cos 2(f> a __ 1 — k 2 sin 2 a/ 2 c os 2 </ > _ A^iy 2 / 1_ ^ 1 


1 — /3 


/x 


" 1 

-fi 


) 

.v 

!_/ 


sin 2 

0\ 

A, 2 '? / 

1 

sin ! 




] — JLi ) 

F 2 / 5 


Therefore 


(z/ 2 — fJ?) ( z^ 2 


1. — ft cos 2(f)' 

^ ) _ AfTf* 




l-^Wl- T i 




- y __ / 2 _ 1 + ff y sin /3 cos /3 cos 7 ‘ (1 — ly) 4 (1 — t, sec 3 /3)= (1 — t, sec 3 7>- v 

\ 1 0/ 


If we write 


<7 = 1(1 + sec 2 fi + sec 2 y) 

II = | (1 + sec 4 y8 + sec 4 y) + \ (sec 2 /3 + sec 2 y + sec 2 /3 sec 2 y), 
this expression, when expanded as far as t, 2 , becomes 


sin /3 cos /I cos 7 




1 

A 8f 2 
A 1 1 1 




II 


The arcs of the three orthogonal curves were denoted dn, dm, df in “ Harmonics,” 
where dn was the outward normal. Since in the present case we are measuring r 
inwards, the element of volume dv must be taken as — dn dm d,f. 

The equations (50) of “ Harmonics ” give 
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do 


cos 6 vdv dO d<j> 


B 


, „ 0 > / 0 1 — /3 cos 2<M /I — /3 cos 2 <t> 

^ t F 

7] ,7J\ 1 + /3 d . a + (3 A* a - /3cos 20G' 

v (y. _ 1 )»(- T —%) cos e ( rz £ - y) I. ! _ p ~ ) 


But 


1 civ ■ 


1-A 

cos 3 /3 cos 3 7 dr 
sin 3 /3 ' A, 3 ]’, 3 


dv _ , s cos /3 cos 7 A ^Fj 3 / 1 

7/A 7 I ”* ^ . * ‘7 /7 „* 0 AH ( T7 ‘ 


i i_ l __ _ \ I _ / _ _ i 1_ n 

dr dO d<j> ~ sin 3 /3 sin 3 7 AT \r x 3 Af*/ L 1 \V ' Fj 3 


1-/3 

, and therefore 

! 1 


Ti — 


a/t 1 ! 3 + G [a^ + iy) “ H )_ 


On comparing this with the expression for p 0 we see that 


dv 

dr 


: p 0 d<T 




V 1 ? 


+ g(-!~ 


W 


I 7 1 O 


1/ 


(5). 


Another form, which will be more generally useful, is found by substituting for r a 
ts value : it is 


dv 


dr dO dcj) 4irp \k { 


3 M ( /d Y [ A-d — I\ 3 ^ cos 3 /3 cos 2 7 


sin 2 7 


sin 2 7 


1 _ r ( 1 _ J 

ov v o? a/, 


4r 


cos 4 /3 cos 4 7 
sin 2 7 


ryA 4 iyA 


1 1 e /1 


r^AY r^y, 


H 




U. 

j Ar 


In order to express this more succinctly let 


<p 


\p = 


6 A t 3 - iy 

7r AT sin 3 7 ’ 

6 cos 3 /3 cos 3 7 


sin 3 7 L r i 4 A i 4 


G 


I? V, 


1_ 
Ar ’ 


^ _ 6 cos 4 /3 cos 4 7 

•——* • Q 

7 T Sill 2 7 


1 _ 1 \ , n( 1 1 

r/vy ryAY /" 1 \ryAy ryy 


// 


( 1 


1 


ryAy r 1 2 A 1 


1 

Ar' 


We note that 


<h. 


l\ da _ 1 M /k 

ddd<f>~ ®7\V 

^ = cos 3 /3cos 3 y h — o' <l5 


Then 


W J 


dv 


dr d6 d(f> 


l-frVC^ - 2 t¥- 4r 3 n] . 

P \' c o/ 

2 L 2 


(O 


( 6 ). 


(«)• 









260 PROFESSOR G. H. DARWIN ON THE STABILITY OF THE PEAR-SHAPED 


The surface r = constant is an ellipsoid similar to J with squares of semi-axes 
reduced in the proportion 1 — 2r to unity. Therefore the volume enclosed between 
the two ellipsoids is 



M (k\* 
P \V 


[1 ~ (1 ~ 2 r )*] 


M 

P 



3 3 

oT 




But taking the limits of 0 and <f> as r,n to 0, so that we integrate through one 
octant and multiply the result by 8, we have another expression for the same thing, 
namely, 

\dv = j (~J | j [<t>r ~ Wr a - f.Qr 3 ] <10 d<f>. 


Therefore equating coefficients of powers of r in the two expressions, 



(7). 


The first of these will be of use hereafter, and all three afford formulae of verification 
in the numerical work. 


§ 4. Determination ofh; Definition of Symbols for Integrals . 

00 o ’ . 

The pear being defined by r = — eS s — 'tfil S’/, with all the f’s of older &, 

excepting f 6 which is zero, we have at the surface of the pear to the fourth order 

r = C- (S 3 f + 2te»f'S a S i ‘ + 

In all the integrations which follow, and especially in the present instance in the 
determination of the volume of the region II, it is important to note that <3>, \P, XI are 
even functions of the angular co-ordinates, and that therefore the integral of any odd 
function of those co-ordinates multiplied by any of these functions will vanish. 
When the odd functions are omitted we may integrate throughout the octant defined 
by the limits |-7r to 0 for 0 and <■/>, and multiply the result by 8. 

Then, only retaining terms as far as e 8 , we may in finding the volume Ii take 

t = — tfdS-R i only even, 

t 3 = e 2 ($ 3 ) fi + 2 'te l fi s S 3 Si, i only odd, 

t 3 = 0. 

To the cubes of small quantities we have, therefore, 
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\ dv = - p J J] [- < PttfSf - V (e 2 (S s f + 2$e>ffS s S?)] dO defy. 

The first term vanishes because Sp is a surface harmonic and <t> dd d<jy is proportional 
to p t) da. 

Thus we are left with 



M / k y 

P \ k oJ 


f f T |> 3 (S s y + 2 te^S,Sf] do defy. 


I now introduce symbols for certain integrals, and in order to bring all the 
definitions together 1 also define several others which will only occur later. 

Let 


c/U' = j| 4> (Sff dO d<fy 
m' = defy 


Pi 


cos 3 /3 cos 3 7 
sin j3 




( 8 ). 


All these integrals vanish unless i is even. For immediate use I also introduce 

ifjp = f j Tosysy dO defy. 


The if/ integrals vanish unless i is odd, but it will appear later that they are not 
actually required. 

I further write 


ff V (S s f dO defy, U = f f O (S,Y dd defy , 


CTa 


6 cos 3 /3 cos 8 7 sin 


sin 3 7 


i^ff fs /_!_ 
JJ L 2 UdiY 


A/iy 


3G 


j._1 \ 

AjTj 4 Adr/jJ 


(M 

AT 


d0d(fy 


( 8 ). 


With this notation we have at once to cubes of small quantities, 
f pdv=-M(~ )V<r, + • ■ 

J r \ / 



But before using this 1 will obtain another integral to the fourth order. It is 
[ rp dv = M (fj J j {i* [e* (S a y + 2Sey?S,S,-> + (2/?S,*)*] 

+f* V W+ *s*ff <W] - «** (*>)*) M rff 


Omitting terms which vanish, amongst which are integrals of the type ^SpS'j, we 
have 
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[ rp dv = {e^ 3 + v {/*)*# + itdfRR - 2e%}. . . (10). 

Returning now to the determination of the mass of + R, and observing that the 
mass of the pear is equal to that of J — R, we have 

M=M[ -J [1 + eV a + 2 tefW]. 

Therefore 

(jf = 1 + eV 2 + 2 tefffl + e 4 S. 

A term e 4 8 of the fourth order has been introduced, hut it will appear that it is 
unnecessary to evaluate it. 

There will be frequent occasion to express /c 5 in terms of kj\ Now 
(jf)“ = 1 - * [«V 8 + 2 tefW + e‘8 - fe‘ (^f]. 

But this will only be needed explicitly as far as e 3 . and to that order 

(*;) = L ~ ••••••••• ( ri )- 


It is, however, necessary to determine f (j 
Now 




k \6 


to the fourth order. 


kV 


I U ~ t [«V 2 + 2 tefW + e*8 ~ |e 4 (cr 2 ) 3 ]} 


S 5 


k \ 5 


fr [eV 3 + 2 + e 4 3 - fed (<r 2 ) 3 ]}. 


Hence to the fourth order 


3. 

2 




9 

1 0 


L ( A 




(Ll). 


It will be observed that the xjj integrals and S have both disappeared. 


5. The Energies %JJ and JR 

If a l9 b l9 Cp are the semi-axes of a Jacobian ellipsoid of mass M l and angular 
velocity co, its lost energy, inclusive of rotation, is 


_JL.1V/J 2 

where SP is the usual auxiliary function. 


'p _j_ b \ + w a 
T 3Jfi 
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The equations to be satisfied by the ellipsoid afford expressions for orb^ and orc^ 
in terms of differentials of AP. If these expressions are added together, or may be 
eliminated, and the expression becomes 

In reverting to the notation adopted here, I remark that |J/, (j^./ will be used to 
denote those functions when the variable is v 0 , and the variable will only be inserted 
explicitly when it has any other value. 

In the present case M r , the mass of the Jacobian ellipsoid, is M and it was 

shown in “ the Pear-shaped Figure ” that 

2 rNf 2 

* = «] ^ = “ J Pi l Qi l 

Hence 

iJ.r = * ¥ £j‘ l9At _ p ^ q . i ]. (12) . 

It was shown in the same paper that the internal potential of the Jacobian 
inclusive of rotation, is 


Wi * + % 


rf'P / x? y% S 2 

^ 77 + &7 7 


Therefore in the present case 


Vj + W(y 2 + * 3 ) 


Mfk \ 3 


2 K\K 


fo<©o — U P^Qi 1 sin 3 /8 (« s sec' 3 y + if sec 3 /3 + 2 s ) 


But the equation to an inequality on the ellipsoid defined by r is in our new notation 
sin 3 /3 (« 3 sec 3 y + y 2 sec 3 /3 + z 2 ) — (1 — 2r); 


therefore 


r, + W (,/ + #) = f f (ff {(ft,®, - P.'QA + 2 rP,>Q 1 >! 

h o \ HJ 


Let us divide this potential into two parts, say U\ TP', of which the first is 
constant and the second a constant multiplied by r. Also let (JR)', (-/A)" be the 
tAvo corresponding portions of the energy JR. 

In order to find (JR)' we have simply to multiply IP by the mass of R considered 
as consisting of positive density. The volume of R is the excess of the volume of J 

above that of the pear; hence the mass of R is M — 1 . 

Therefore 
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Subtracting this from \JJ as given in (12), 


iJJ - (JR)' = 


if 3 


(fo®o 



But the latter factor was found in (11) as equal to fo — |-e 4 (erf 1 . The term fo 
only contributes a constant to the whole energy and may therefore be dropped. 
Accordingly 

= . . . os). 

/v 0 


For the other portion (JR)" we have 



Then by means of (10) 

(JR)" = f U"pdv = 3 f (~Y P/Q, 1 f r P dv 
J r li {) \h\f J r 

/If 2 / Z» \5 

= I I- r P/Qi 1 + >- (./?)-'/>;’ + - 2e 4 y (14). 

# () \A 0 / 


In the terms of the fourth order we may put (k/kf equal to unity. Therefore 
combining (13) and (14) 


w-™= -tf(^ 


J/ 3 


^P^QxVs + l y (f 0 ® 0 - Pi 1 Q 1 1 )(«r a )» + 

A '0 

-4P,>Q,'S«*/ j w - P^Q/S (/•«.') (15). 


6. Surface Density of Concentration C; Energy CR. 


The region R being filled with positive volume density p, is concentrated along 
orthogonal tubes on to J, and there gives surface density S. 

To the first order, by (5), 


dv 

ch 


peeler 


Integrating with respect to 
small quantities 


2r 


cos 3 /3 cos 3 7 / 1 




Ay 1 F 3 


AfTy \iiy i y 

t from the pear to .7, we have as far as squares of 


8 = ~ Pop 


eS s + %f s Si + e 


2 cos 3 /3 cos 3 7 


A.Ty \\ 


h + 7 - c>)(s,f 
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It is now necessary to express S /p 0 in surface harmonics. The first two terms are 
already in the required form; for the remainder let 


trji’Si* 


cos 3 /3 cos 3 7 /I 1 


Ti + ri-0 KW 


A^Fj 3 \Aj 3 1 l 1 , 3 

Multiplying both sides by Sj s ®clOd<f> and integrating, we have 

’lity = COS 2 £ cos 2 y + l !, - g) (S & )~S; s cW d(j) 

= fj ' v (s s ySi* dd d<f> = co;\ 

Therefore rji * = */<£**, and vanishes unless i is even. 

When i = 0, rj 0 = and since bj (7) <j> 0 = 3 , and co 0 = 'P (S s )*d0d<f> 

we have r} 0 = . 

Hence we have 

= ~~ Pop e $s + 3 e2(T 2 + | ( e2 ^ $ 

This is expressed in surface harmonics, the middle term being of order zero. 

By (51) of “ Harmonics” the internal potential of 8 is 

<#3 (0 ©3 (-0) *3 + - b^Mo + 2 U ~ +//) p/ (v) ©/ (v 0 ) # l 


<T 


2> 


7.= 


4 \ h / 

We have P/(^) = 




cos 3 /3 cos 3 7 1 


^_ T g ut n e f ore 

2v 0 clv 0 sin/3 A, 3 I \ 3 ' 


proceeding to use this I will introduce a new abridgment, and let 


% s = f m s 


b =<&/ 


dVn 


(16). 


Then 
V = - 3 


U (») ® (v„) = &• - r ”, and 


Mlk_ 
lc 0 \ . k 0 


e ^3^3 T" b 2(r : a^o + ^ ( e 2 — + //) 'U//?/ 


+ 3 


AT / y cos 3 /3 cos 3 7 


*0 \*( 


sin /3 


& 


eB 3 + 2 (e a ^ + j?) Ark 




j A/T\ s 


In order to find the energy CR we multiply F„ by the element of mass 

p dv —■ \M ( —) [cl> — 2rT r ] dr d9 d<fi, 

\ ,c o/ 


and integrate throughout R. 
VOL, cc.—A. 


2 M 
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Now 


V c p dv 


3 _ 

8 h \ b 
/Tn \ A'/ 


^ 4“ 3^ 2(r ; 2^0® + S (^ 2 — + jO j 


i- 


0 Ml /AY T ] e^3^ s .+ ieV,^F + s ( e 2 g + ft ) a,**# 

<Pi 


h \V 


jlAYAY T cos ; 


+ 8-T- l IT ! T 
/; 0 \ /b 0 ' 


AtAA La* _A_ + v L ^ . 

sin £ 1 “ 3 A 1 3 r i 3 h Z \ 4>f + ‘ 




Sf 


Let us integrate these three lines separately. 
First integral 

m 1 k V? rf \e® s ®S 3 + + t U 





J { + (A ) 3 


w 


1 VU 


?)<&•*&} {eS, + t/m <16(14 


Second integral 
M 3 / k \B 



e® 3 VS 3 + le*<r 3 ®oV + t ( 4 ^ + //) 


{^(S.Y + 2 $efiS s S?} cie d(f> 


M* fk_ Y 
*o \ *o/ 


(<r 2 ) s + 2 |.e* +>./?«,* j &* + 2^^/ 


Third integral 

df 2 / 1c \ r> cos 2 /3 cos 3 7 C j" 
sin /3 


2 *o \*b/ 


Al_ ^ As _l /’.A as»—A 


* 3 Adr-d^^V ^a^j 

(e 3 (4 ; ) 2 + 22 e/^S/} <16 <14 


M2 ( b \5 f / \ 

= f~(f- LK?. + e A P/2S/ + 224/?B 3 L 

Y) \ /v o / l \ 9* / 

All the terms, excepting the first of the first integral, are of the fourth order, and 
in them we may put (&//%) 5 equal to unity. 

Therefore 

if 2 / k \ 5 

CR = 3~(~) e%</> 3 
lv 0 \ /l/ 0 / 

+1 {4 r§a 0 ( 0,) 3 + 2 (2a>/+ wpt) 

+ tej ![4 (%• + 0 ^) <»; + (w + m 3 ) P n + 22 : (fifm 4 i\. 


(17). 
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§ 7. The Energy \GG; Result for \JJ - JR + CR - \GG. 

From the last section it appears that the potential of G at the surface, where 


0, is 


F n 


g M/kV 


hVv. 


+ ieVA + t (e 3 + //) <&S t ‘ 


For the mass of an element of the surface density we have 


'k\ 3 


WO/ 


p 0 8dcr =-iM{ J cp \ eS 3 + ieV 3 + S (e 2 ^ + //) S/ f # 
* / ' 


0/ 


These are to be multiplied together and half the product is to be integrated. 
Then bearing in mind that | j <E> d6 d<j> = 3, we have 




+ l e4 K)% + S ^ + fi 


In the terms of the fourth order we put (Jc/k Q ) n equal to unity; thus 


i CO=s s mv 

2 2 *oU. 


<*aw» + fr{« 


l L 




( w/) 2 

J 


+ 2te?f s a>M i s -f tifryvity . 


Combining this with (17) 


OR - ICC: 


= f f o (fJ « 2 a 3 c/>3 +1 ^ {e 1 [Wo fo) 3 + s (av + 3s ; >/) 


ft)/ 


</>/ 


+ te*ff2(®: + 2® 3 )a>° + (w + 2® 3 ) P fi + $(f°ymn ■ (is) 


We are in a position to collect together all the results obtained up to this point. 
Now \ JJ —JR, as given in (15), contains Pi'Q, 1 , H ie latter of these is what 

is now written a> and since the ellipsoid is critical F 5 1 1 Q 1 1 = = a* 

Collecting terms we find that the terms of the second order disappear, and that 


w 


JR+CR- \GG= | , 

k 


M 2 


; (i (cr 2 ) 3 + 2 Q +1 (av + Wpi) 


H 


+ [2av + (is/ + 2b 3 ) pf\ -1 (f s f (a - a s ) h ) • (19). 


The reader will recognise that the last term involves the coefficient of stability for 
the deformation S'/. It is important to note that if S’/ is of odd order there is no 
term with coefficient e % f\ 


2 M 2 
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8. The Term ~~ \MiVcT. 


In the Jacobian ellipsoid 


bf + c 


1 2 
ar 


v + !«i 


(W 

da. 


In the present notation this is 
(1 + cos 3 /3\ 2 

Hence 


3 M \ sin 3 j3 j h 


(Sofio - IPi'Qi') 




3 


(a, - fa 3 ). 

(Mdf 


1 + cos ! 

I now make the following definition 

Si = sin 6 (1 — k ' 2 cos 3 (f> )% 


m 


¥' 


so that 
Then 


z = kS v 


Md 


— zp dv = zp dv — zp dv — — zp 
J ?_#* J j Jr J r 


= - M (~Jk [ ff [<1> - 2 r'P] ^ dr d6 d<j> 

= i¥ /|Y* (f [«£ (e/S, + 2/Td) + 'Pe 3 (^ 3 ) a ] ,S\ dd d<j> 

= M (£) S kM, 

\ “'A / 


Therefore to the required order 


= 


3 M 3 sin 3 /3 , 

l”+ cos 3 /3 




( 20 ). 


We again note that this term in the energy does not introduce any term with a 
coefficient e~J]. Hence thus far the whole energy for harmonic deformations of odd 
order is of the form Le 4 + M (fd) 2 . 


§ 9. Double Layers . 

It remains to determine the value of \DD in the energy, and for this purpose we 
must consider double layers, according to the ingenious method devised by 
M. Poincare. 

Let a closed surface S be intersected at every point by a member of a family of 
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curves, and let a be the angle between the curve and the outward normal at any 
point. At every point of S measure along the curve an infinitesimal arc r, and let r 
be a function of the two co-ordinates which determine position on S. The extremities 
of these arcs define a second surface S', and every element of area da of S has its 
corresponding element da' on S'. Suppose that S is coated with surface density 8, and 
that S' is coated with surface density — S', where 8 da = S' da'. The system SB' 
may then be called a double layer, and its total mass is zero. We are to discuss the 
potential of such a system. 

Let U (~\-) and U ( —) be the external and internal potentials of density 8 on S, 
and U 0 their common value at a point P of S’. At P take a system of rectangular 
axes, n being along the outward normal, and s and t mutually at right angles in the 
tangent plane. 

In the neighbourhood of P 


U(+) 


u » + 45(+) + 4;(+) + 4t(+) 


ds 


dt 


U(- 




ds 


dt 


In the first of these n is necessarily positive, in the second negative. 

Now ( + ) = ~ (—) = ; and the like holds for the differentials with respect 

ds v ‘ ds v ' ds 1 

to t. 

Also by Poisson’s equation 


dU, . dU , . s 

(“) “ ife( + ) = 4,rS - 


dn 


Let PP r be one of the family of curves whereby the double layer is defined, and 
let P f lie on S\ so that PP' is r. By the definition of a the normal elevation of S ; 
above S is r cos a. 

Let v, v' be the potentials of the double layer at P and at P\ 

The potential of S f at P f differs infinitely little in magnitude, but is of the opposite 
sign from that of S at P ; it is therefore — U 0 . The point P f lies on the positive 
side of S at a point whose co-ordinates may be taken to be 

n = r cos a, 5 = r sin a, t = 0. 

Therefore the potential of S at P f is 

rr , dUt , x . . dU 

' 7 » +Tcos “«fa( + )+ 


/_ d U , - v | 

v = r cos a —“ (+) + r sin 
dn v 1 


a 


dU 

ds 


Therefore 
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Again the potential of S at P is Z7 0 , and since P lies on the negative side of S' and 
has co-ordinates relatively to the n, s, t axes at P' given by 


n = t cos a, s = r sm a, 
since further the density on S' is negative, we have 


Therefore 


t 0 


v — a = r cos a 


dn 


dU, 

dn ' 

\ . * (W 
) + r sm a — . 

7 ds 

T 

- 7 -( + ) = forS cos a. 
dn x ’ 


The differential with respect to n of the potential of S falls abruptly by 47rS as we 
cross S normally from the negative to the positive side; and the differential of the 
potential of S' rises abruptly by the same amount as we pass on across S'. It 
follows that clv/dn on the inside of S is continuous with its value on the outside of S'. 

The surface S to which this theorem is to be applied is a slightly deformed 
ellipsoid, and the curves are the intersection of the two quadrics confocal with the 
ellipsoid which is deformed. The curves start normally to the ellipsoid, and where 
they meet S the angle a will be proportional to the deformation whereby S is derived 
from the ellipsoid. It follows that cos a will only differ from unity by a term 
proportional to the square of the deformation, and as it is only necessary to retain 
terms of the order of the first power of the deformation, we may treat cos a as unity. 

We thus have the result 

v — v f = 


Suppose the curve PP' produced both ways, and that M 0i M l are two points on it 
either both on the same side or on opposite sides of the double layer. 

Let M () M 1 be equal to £, let £ be measured in the same direction as a, and let £ be 
a small quantity whose first power is to be retained in the results. 

Let v {)i t\ be the potential of the double layer at M 0 and M x respectively. 

When £ does not cut the layer we have 


y dv 

”o - = - ( dn • 


and when it does cut the layer 


dv 

v 0 - V] ^4 rnb-i-. 


In the application which I shall make of this result the surface S' will actually be 
inside S. Then v 0 will denote the potential at any point not lying in the infinitely 
small space between S and S\ and v 1 is the potential at a point more towards the 
inside of the ellipsoid by a distance £ ; 8 is the surface density on the external surface 
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S and 7 is measured inwards. If then we still choose to measure v outwards, as I 
shall do, our formula becomes 

% ~ v i - or 0, 

according as £ does or does not cut the double layer. 

It may be well to remark that v being proportional to r8, £ dr/dn is small compared 
with 4 tttS. It is also important to notice that the term 47 ttS is independent of the 
form of the surface, and that dv/dn will be the same to the first order of small 
quantities for a slightly^ deformed ellipsoid as for the ellipsoid itself. 

We have now to apply these results to our problem. 

The position of a point in the region R may be defined by the distance measured 
inwards from J along one of the curves orthogonal to J. The surface of the pear as 
defined in this way is given by e, a function of 6 and <fi. Any point on a curve may 
then be defined by se, where s is a proper fraction. If s is the same at every point 
the surface s is a deformed ellipsoid ; s = 1 gives the pear and s = 0 the ellipsoid J. 

If dcr is an element of area of J, the corresponding element on the surface s will be 
(1 ■— Xes) dcr. The value of X will be determined hereafter, and it is only necessary 
to remark that it is positive because the areas must decrease as we travel inwards. 

Let s and s + ds be two adjacent surfaces ; then the mass of negative density 
enclosed between them in the tube of which (1 — Xes) dcr and (1 — Xe {s + ds)) dcr 
are the ends is — pe (1 — Xes) clcr ds. If this element of mass be regarded as 
surface density on s, that surface density is clearly — peels. If the same element of 
mass were carried along the orthogonal tube and deposited as surface density on J, 
that surface density would be — pe (1 — Xes). The sum for all values of s of all such 
transportals would constitute the condensation — C already considered. 

The double system D consists of the volume density *— p in R y and the positive 
condensation -j~ C on J, the total mass being zero. 

Let z, a proper fraction, define a surface between J and the pear. Consider one of 
the orthogonal curves, and let V {) be the potential of D at the point P where the 
curve leaves J and V z the potential at the point Q where it cuts z. Then I require 
to find V 0 — V z . 

Since s denotes a surface intermediate between J and the pear, -- ( V 0 — Vd) ds is 

the excess of the potential at P above that at Q of surface density — pe ds on s and 
surface density + pc (1 — Xes) ds on J. Such a system is a double layer, but there is 

a finite distance between the two surfaces, and the form of ~ ( V 0 — V g ) will clearly 

be different according as z is greater or less than s. 

The arc es may be equally divided by a large number of surfaces, and we may take 
t to define any one of them. Now we may clothe each intermediate surface t with 
equal and opposite surface densities zlz pz [1 —• Xe (s — t)~\ dt> 
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The density + pe [1 — Ae (s — £)] dt on t , together with — pe [1 — Ae (s — t — dt)~] dt 
on t + dt> constitute an infinitesimal double layer; and since the positive density on 
each t surface may be coupled with the negative density on the next interior surface, 
the finite double layer may be built up from a number of infinitesimal double layers. 

Hence (F 0 — Vi) dt dt is the excess of the potential at P above that at Q of an 

infinitesimal double layer of thickness edt , and with surface density pe [1 — Ae (s — 1)~] dt 
on its exterior surface. 

We may now apply the result v Q — v l — = TttSt or 0, according as £ does or 

Cl/lb 

does not cut the double layer, and it is clear that 


cP 

ds dt 




= 4?rp€ 3 [1 


— Ae (s — £)] or 0, 


according as z is greater or less than t. 

In the next place, we must integrate this from t = s to t = 0, and the result will 
have two forms. 

First , suppose z > s; then for all the values of t, z > t 5 and the first alternative 
holds good. Therefore 


d 

ds 



V — 

* z 



4? rpe 2 [»S' — -|AeF] . 


Secondly , suppose z < s; then from t = s to t = z, z < t and the second alternative 
holds, while from t = z to t = 0,z>t and the first holds. Therefore 

l( F »- V ’~ 62 d £) = ^ 0 s - - i 2 ’)] • 

We have now to integrate again from s = 1 to s = 0. 

From 5 = 1 to s = z 9 z < s and the second form is applicable; from s = z to 5 = 0, 
z > s and the first form applies. 

Therefore 

dv c 1 r* 

V 0 — V z — ez ^ = 4.7Tpe 2 j [z — \e (sz — |-£ 3 )] c/s + 47rpe 3 j [5 — JAes 2 ] ris 

= irrpe 2 {z(l — z) —■ Ae [^2 (1 — 2 3 ) — |rZ 3 (1 — 2 )] + -|2 3 — ^-Ae2 8 } 
= 2npe 2 {2z ~~ z 2 ~~ Ae (2 — 2 3 + ^ 3 )}. 

Finally, we have to multiply — 4 ( F 0 — V g ) by an element of negative mass at the 
point defined by 2 and integrate throughout II. The physical meaning of this 
integral will be considered subsequently. 

We have already seen that such an element of mass is given by 
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•—pdv == —- pc (1 — Xez) da dz 


and the limits of integration are z = 1 to z = 0. 
Therefore 




i dv 


rr ff dV 

j e 3 (1 -- Xez) { 2 z — — Xc ( 2 ; — £ 2 + } dz da + i/> | € 2 2 (1 — Xe^) — d£ dcr. 


: 7T/T 


c?F . 


In this expression we neglect terms of the order e 5 and note that e 3 £ 2 —- is ot 


that order. 
Thus 


(F 0 — F,)pda^7rp 2 ffe^ dzda (s= 1 to 0), 

= §7/yr fe 3 (l — Xe) da + ip je 2 “ da, 


the integrals being taken all over the surface of the ellipsoid. 

We must now consider the meaning of the integral -J | (F 0 — F,) p d?\ 

Let P be a point on J and Q a point in It on the same orthogonal curve. 

Let — U be the potential at Q of the density — p throughout 11, and — U 0 its 
value at P. 


Let 8 be the surface density of the positive concentration on J, W its potential at 
Q, and IF 0 its value at P. 

The lost energy of the double system consisting of — p throughout R ? and 8 
on J is 


if Updv + i 


W 0 8 dcr 




U 0 8 da 


Wp dv. 


This is equal to 

i f ( U - W) p dv - j- f( U 0 - W 0 ) 8 da. 

Consider the triple integral (CT 0 — W 0 ) pdv. Here dv = e (1 — kcs) da ds ; also 
U 0 — W 0 is not a function of s, and the limits of s are 1 to zero. Therefore 


( U 0 — W 0 ) pdv = ( U 0 — TF 0 ) e (I — Xes) p ds 

J J J J J J o 


dcr. 


But e (1 — Xes) p ds is equal to 8 the surface density of concentration. Therefore 
Jo 

f f[»» - K] 8 do- = f j|(£A 0 - W a ) P dv. 

VOL. CC.—A. 2 N 




274 PROFESSOR G. H. DARWIN ON THE STABILITY OF THE PEAR-SHAPED 

We may now revert to the Gaussian notation with single integral sign, and we 
see that the lost energy of the system is 

ifKB'o- U { )-(W-U)-\pdv. 

But W — U is the potential of the double system at Q, and is therefore V e ; and 
Wo - U 0 is the potential of the double system at P, and is therefore V 0 . 

Accordingly the lost energy 

iDD = i\(V 0 -V s )pdv 

.( 2i ). 


10 . Determination of e and X. 


e is the arc of the orthogonal curve from J to the pear. 

The arc of outward normal is connected with p and our variable r by the equation 


dn : 


¥ 

P 


vdv — — dr. 
P 


it follows that 


p 0 ~ 7 ° dr, integrated from J to the pear. 

j p 


By (50) of “ Harmonics,” with the notation of § 3 of this paper 


k 


A. \ i - / 


sin y 3 


(1 .iWi- 

- - T iV 

\ Ad/ 1 

JV/ 


1 + 0 


P V ( V 2 — 1)* ( v~ 

K \ 1-/3 

Therefore 


cos y 3 cos 7 (1 — Tj)* (1 — r x sec 2 / 3 )* (1 — sec 2 7)* 


Vo 

P 


Ad 


1 - - J 

W 


(1 — Tp)2 (1 — sec” f 3 f (1 — r x sec 3 7)* 


2*1 


— 4 . - - 

A x 3 ‘ If 


2 G\ 


= 1 —- T 


cos 3 /3 cos 3 7/I 


Ad iy 


Adi? 

Integrating this from J to the pear 

e = — jPo [b>3 + + i e 3 - 


2 6r ). 


cos 3 /3 cos 3 7/I , 1 

~~A/f/ \V + f7 3 


)(«„)’ 


. (22). 


We have, moreover, by the formula before integration 
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cln = p 0 


_ 12 s ! Z 3 cos2 7 /J- , JL __ on 

T A , 2 !’, 2 \Ai 2 ‘ Tj 3 


dr . 


Also to the order, zero — n = p^r. 

Since — n is what was denoted in § 9 by e«s*, the element of volume is 
- (l + Xn) dndcr , and this is equal to 


P o [1 - \Po T ] 


or 


Ih 


i 


cos 2 /3 cos 2 7 M 1 n \l J , 

T AjV ( a l 2 + ly 2 /J ^ T ’ 

cos 3 /3 cos 3 7/1 1 


1 T { + A, 2 ! 1 ! 2 ’ hi 2 i\* 2G 


dor clr. 


I 1 cos2 ^ cosii V ('JL _j_ A _ 2G 


But by (5) the element of volume is 

' p<) [/ a/iv VAj 2 1 17 
Equating coefficients of r in the two expressions we find 


_ cos 3 )6 cos 3 7/II. 

Pu&il-Y Uf 1?'* ‘ * 


(23). 


§ 11. The Energy §7r/> 3 j e 3 (1 — \e) dcr. 

From (22) and (23) we have 

ft 3 1> w + (S,y S> + |e‘ Q, + A 




2G)(s a y 


So that 

e 3 (l — Xe) = - p 0 * 
Again from (6) 

=i^p 
Therefore 
1 7 rp 3 e 3 (1 — Xe) dcr 


<W + «*VW«*+«*11(A + h)- 3s»«» 


A/Tj 2 1 2 \A X 2 n 17 


2 _ 2 , Af / ft Y cos 2 /3 cos 2 7 (16 del) , if 2 / ft Y 

‘.“ <lv 7 » - .• • 2 ~ 1 (i /. I t, ) 

/l 'l) \ A, o/ 


tl> 


8 p v*< 


H>/ 


sin 2 /3 AfTj 3 


cos /3 cos y sin ft — </# d<j>. 

Ayl l 


1 ^V-Y 

16 7 , It,/ 

/l/ o W 


cos cos y sin /3 


, 8Sa ,,.*W# 

L a/t/ ^ Ji Afi'y 




Ad'IY 

2 n 2 


1 . 1 - 


. A 2 1 P! 

\“l 1 1 


3(7 


d9 d6. 
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When this is integrated we may put (k/k v f equal to unity. In the integral the 

3 If 2 sin 2 /'I 


first term vanishes, and the second term gives 
third term we substitute for 4> its value and have 


2 k 0 cos /3 cos 7 


%e l f s p-\ In the 


3 M 2 cos 3 ft cos 3 7 sin ft 4 ff 1 /l 1 \ f 5 / 1 1 \ 0 H dOdrj> 

” it k 0 .sin 3 7 C J j A/rd ~ V/I 2 \I? + A?/ ~~ J AT ’ 

which is equal to 

1 If 3 ^ 6 cos 3 /3 cos 3 7 sin ft f f [ 5 / 1 1 \ 0 ,- Y / 1 1 \ ] defy 

tv sinDy JJ j * \f?V “ IVA/V "™ ' itT vF-fAf 2 ” J AF 


By the definition ( 8 ) this is equal to — — e 4 , cr 4 , 

/u 0 

Hence the required term in the energy is 


3 M 3 

2 £ 
/L 0 


cos ft cos 7 




(24). 


§ 12 . 


The Energy |p 


o rfF , 
an 


It is first necessary to determine dVIdn . 

Suppose that the ellipsoid J is coated with surface density 8 , and that a second 
surface is drawn inside J at an infinitesimal distance t, and coated with negative 
surface density — S ; ? so that the two form a double layer. Then rS being a function 
of the two angular co-ordinates on the ellipsoid may be expanded in surface 
harmonics ; suppose then that 

rS^th'S*. 

o 

Consider the two functions 


V e = % 4tt/C (V * — I ) 1 ( 7 ' 0 2 


Vi 


1 + pf dffi(v 0 ) 
1 — ft dv 0 


(p) for external space, 


f}/ (p) for internal 

(iVa 


space. 


Since these functions are solid harmonics, the matter of which V e and Vi are the 
potentials is entirely confined to the surface of the ellipsoid, and since they are not 
continuous with one another, the ellipsoid must be a double layer. 


Now 




dv 


PP/MTV - i) J 7 


t+ 

:i - '$} 


and therefore 
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/ 1 V i / 2 1 + 

("o-i >\ v °-r+p. 


Hence at the surface of the ellipsoid 

V e - Vi = S 4 ttA/£/ = 4 tttS. 


But this is the law found in § 9 for the change of potential in crossing a double 
layer, and hence F 0 Vi are the external and internal potentials of the double layer 


tS. 

Since 


d 20 d 

dn W vdv 5 


C 1L 

dn 


dV A 

dn 


d,V « 4ny 0 , 2 
“7” = 2 (V 

dn o 



t + £\* , , <£|l/ d@/ ~ 
1 - / 3 / * ' dz / 0 * 


(25). 


This result will hold good to the first order of small quantities if the surface be a 
slightly deformed ellipsoid, such as was the surface defined by ^ in § 9. 

In the elementary double layer t the density was pe [1 — \e (s —* t)~\ dt , and 
the thickness was e dt, so that the thickness multiplied by the density was 
pe 2 [1 —• \e (s — t)~] dt dt . Since, however, we only need this to the first order, we 
may take it as pe 2 dt dt. It will now be convenient to change the meaning of hf to 
some extent, and to write 

e 2 = 2 lx* Sf. 

o 

Thus for the elementary double layer we have 

r 8 ~ pdtdt% hi s . 

0 


d V 

It follows that in applying the formula (25) to determine for the double 

CLTb 

system D, we may say that 

dadl dn f *■ 0 ' V° 1 -/ 3 / ' * 0 " 

Since the right-hand side does not contain t, we have only to consider the integral 

O ds dt = f s ds = ly. 

0 Jo 


Thus, for the system D, 


dV_ 

dn 





1 + / 8 \* 

1-/3/ dv 0 dv 0 


( 26 ). 





278 PROFESSOR G, H. DARWIN ON THE STABILITY OF THE PEAR-SHAPED 


This result may also be obtained as follows :—-To the first order we may concentrate 
the negative density in the region R on a surface bisecting that region. We may 
then consider the positive concentration C on J, and the negative concentration on 
the bisecting surface as an infinitesimal double layer of thickness -|e. We have seen 
that the surface density + C is — ppe$& and that € — — peS < 6 (in both cases to the 
first order only). Thus the density S of + C is pe, and the thickness r of our layer 

is ^e; the product therefore rS is IpeL 

00 

Hence rS = ±pe 2 = and thus we arrive at the same result as before. 

0 

I now introduce an abridged notation analogous to that used previously, and write 


w = 


dv 0 clv {) 


We then have by ( 26 ) on the last page 


^ = s , 

an o p 


where 




By (22) to the first order 


Assume then 


°~'D 2 (S B — e'Vc 2 c V eos dV. 
AA IV) sin2/3 Al AV 


cos 3 p cos 3 y (S s f " s „ s 
sin p A x 3 f o “ 1 ' 


Multiplying by $>S* and integrating, we have 


pi = 


Hence 


00 n 9 p?,l3 n * 

c 8 = — %~Si, and therefore hi — - 7 . 

sin 0 o </>/ sin p <f>f 


Substituting in ( 26 ) 


c os p cos 7 pj , t v . 


dn “ ^ 2wpope “ sin® /S 


2 | - e*p Q sin P Y 20/$*. 

0 # (pi 


plS CO Q ? 

c 2 z=: ——- % ~~ Si 9 and 
Sill 0 0 <fri 
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Since on integration the terms involving products of unlike harmonics will dis¬ 
appear, we have, as far as material, 


clV 

dn 


_ 3 . 


M 




5 \2 


Now 


PI 

0 V&\ 


kV 


3®./ (s?y. 


ipPo d<r = ~h M (pr) de 


Since the term which is being determined is of the fourth order in e, we may put 
k/k 0 = 1, and we have 


ip f«* f *=* f ,'*! ff( |)' ** w 

M * . s (^) 2 


a --- e 4 X 


■mi 


(27). 


*0 0 4 >i 

Since “ 1 ? IPo ( v ) = 0 and 39 0 = 0, the term in S corresponding to i = 0 

vanishes. 


§ 13 . Terms in the Energy Depending on the Moment of Inertia . 

We have to determine A n the moment of inertia of the region R considered, as 
filled with positive density. 

In order to obtain this result, we must express if +$ hi terms of surface 
harmonics. This was done in § 12 of “Harmonics,” but as a different definition of 
S 2 and S% was adopted there from that which I shall use here, it is easier to proceed 
ah initio . 

Let D 2 = 1 - and 

fa>) ! = 10 +W = «i + * s + D). 

For both the suffixes 0 and 2, we have (f -f- </ 2 = 1, and 


. 2 - 3 r/ 


k ' 2 = q ’ 2 


1 . — ■)(£' 

1 _ f ... Yc? ’ 7 ” 7 “ 1 - 2 (f 

In accordance with equation (10) of “The Pear-shaped Figure” I define the 
harmonics as follows :— 

S 2 = (V sin 2 0 — q 0 2 ) OV — K ' 2 cos 2 <£), 

S / — (V sin 2 0 — g 2 2 ) (g' 2 2 — k ' 2 cos 2 <£)• 

Now y~ = /c 3 (r 3 — 1) cos 3 6* sin 3 <p, z 3 = /cV 2 sin 3 d (1 — /<' 3 cos 3 </>), 

and, ” 3 


k 2 - <f = g' 3 - k' 2 


{ t<P 


1 — t, 2t cos 3 /9 cos 2 7 


v~ — — , 1 , where r x 


sin 3 /3 : 


Aprp 
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— cos 2 /3 — - Tj + si.ir /3 sin 2 0 — (cos 2 /3 — r,) cos 2 c/> -|- (cos 2 /3 — A 2 — kVj) sin 3 0 cos 2 <£. 
Let us assume, as we know to be justifiable, 

^(^ + 3 *) = A^ + RS^ + C 

- - [A 2u Yo 3 + BYW - CTJ + [Ar/o 2 + B3?] K 2 sin 2 0 
+ [Ag'o 2 + Bg 2 2 ] k' 2 cos 2 <j> - [A + B] kV 2 sin 2 0 cos 2 <f>. 

If we equate the coefficients of sin 2 0 and cos 3 <jf> in these two expressions, we have 


a s v + b ? '/=^ 


a j_ TL» a — — c -°- 2 |S 


A? 0 8 + .B? S 8 


The solution of these equations may be written 


1 A cow 2 /3 


D + *»/ 2£g' 0 s 


B = — — 1 - /1 + 

22? 2s * V r d - **/ ^ 2/yt/ 


The simplest way of finding C is to put sin 2 6 


cos 2 </> = Ya > 80 that 

fC 


$ 2 3 = 0 ; we thus find 


Now for brevity write 


C = i(l + COS*0) - frj 


4D £ 0 3 cos /? cos 7 \ 1) - \ 


cos 3 /3 


cos /3 cos 7 \ " 1 1) 


cos 8 /3 


We then have 


» ^ cos 8 cos/y r 

A = 2 -.W;r A 




B = - 2 


cos /3 cos 7 




Hence, substituting for 77 its value, 


y 2 -f $ 0 cos /3 cos 7 


2 — W 7 

sir L 


1 + COS 3 8 


fl ( LB,-Man+ ™ 


COS 2 y3 cos 2 7/1 N 3 1 N/ 4 1 ) ' 

' D sin 2 /3 W A/I? AfTi ’* 


_ 1 1/ [l) S (<p _ 2 tT') 

dTded4~ sM \kj 1 


Therefore 
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(if + 2 2 ) p civ 
dr (16 d<j> 




iY 12 cos ?-™P(L<S->S 2 - M<t>Sf) + ® 


sin 3 /S 


4r 


cos /3 cos 7 
sin 3 /3 


3 sin 3 j3 

{IWS 2 - M9SJ) - 2r 


1 -f- cos 3 /3 

-L_ 

3 sin 3 /3 


cos 3 /3 cos 3 7/I 1 <E>$ 2 


i)shF/3 \ 2 WIV ffVA^f 


h + 


M- 


When we integrate throughout the region R the limits of r are — eS a — "tffSi to 


zero. 


Accordingly 


A r 


Mkf 


0 ‘ff 


cos p cor 7 _ }/(j) Q 2\ _j_ 1 <p 

sin 3 /3 1 2 3 ’ ■ 3 sin 3 0 


eS, + $f?S° 

2 


+ 


_ COS /3 COS 7 , T T ~ - r c e>\ . 1 + COS 2 P 

2-r™ (LTlSo — AfToSi, 2 ) + —-r-~ T 


sin 8 /3 


3 sin 3 /3 


e 3 (S') 2 


+ 


■ l/V 


cos 3 /3 cos 3 7 
2j 9 sin 3 /3 


1 <M, 1 CRN/ A <l> 

2T1 2 ~T 3 V A 


/o 2 a/w 2 v AjTj 2 T 3 - a^f, 2 






+ 


_/ P%. __ pi 1 4 


2D sin /3 \q f ^ (// 


h + 3 d Po ) ^. 


The moment of inertia of the ellipsoid J is 


Aj = }M( YY/0 2 1 + c f ^ = Alkf 


sin 3 /3 


‘1 + cos 3 ft ^ 1 + eos 3 /3 


5 sin 3 /3 


3 sin 3 /3 


CTo 


Also 


MW. 


,¥ 3 3 sin 3 /3 3if 3 1 sin 3 /3 


A * 4;7r p cos yS cos 7 2£ 0 * 2tt p * cos , 8 cos 7 ’ 


Lastly, to the required order we may put (/c//c 0 ) 5 equal to unity in the expression 
for A r . 

Then 


i(Aj-A r ) §o>* 


3if 3 So) 3 [(1 + cos 3 /3) sin. /3 


2k 0 2nrp [ 10 cos /3 cos 7 


- + L (Ah + *<h) - if WW + e V) 


+ 


__ 3 __ s ^ n2 & (£2 _ A 2 


4D cos /3 cos 7 \ i/ 0 2 


h + i-D* 


( 28 ). 


This completes the expression for the lost energy E of the system, which may now 
be collected from ( 19 ), ( 20 ), ( 24 ), ( 27 ), and ( 28 ). 

VOL. CC.— A. 2 O 
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§14. The Lost Energy of the System; Solution of the Problem. 

If the several contributions to the energy be examined, it will be seen that if i, the 
order of harmonics in fiSi, is odd, there is no term with coefficient e 9 f? in E; this 
follows from the fact that the co and p integrals vanish for the odd harmonics. Hence, 
as far as concerns the odd harmonics, E involves f s only in the form {fff. The 
condition that the pear shall be a level surface is that E shall be stationary for 
variations of the f’s and of e. It follows that when i is odd f s is zero. We may 
therefore drop all the odd harmonics, inclusive of f, and it is clear that the term 
— \M<P. a) 2 in E (given in (20)) vanishes to our order of approximation. 

For the sake of brevity, I adopt a single symbol for the coefficients of the several 
kinds of terms in E Therefore let 


d. 0 — % A [j- (cr 2 ) a + 2£J + | + I h s pi s ) f s ~ 3 _cr 4. T" 4 ^ ' 

sin ' 15 /3 


2Bi “ ~ cos/S cos>y P*’ 

(1 + cos 3 ft) sin ft 


10 cos ft cos 7 J 

& = 1 / 0 )% — M(i)% + 
t — L<j > 2 , 


sin 2 yS 


Pi PJi I 4 7) V 

f 2 + 3 JJ P0j > 


2 " r 4 D cos ft cos 7 \q f Q q'% 3 


= where %■ = •PtW, 38 / 

With this notation 


©/ 


"0 


dv 0 cIvq 


E 


t {** +*f Byf! - ? w +!”p < a +^+ & -w)}- 


Let us now make E stationary for variations of e and f s . 
First, by the variation of any f? excepting f 2 and we have 


13 s 

St = Le 3 


0 ? 


. . (29). 


On eliminating all these fi, we have 


E _ 3 


A„ + Z ®) «* + 2B a e% + UWff - C, (/,)* - Of (/,*)* 


+ \ 2 “ (« + k* + e/, — . 
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By the variations of / s , and e 3 , we have 

B+ ~ t = 0, B,V - Cfff -~t = 0, 

( A o + i ( -jx) e ' ! + Bif s , + B,%<‘ + = 

But from the first two of these equations 


Therefore 


j> f _ (A) a 2 i ! A c 2 y 

-«2/3- Cf e ' 4 ttp C' 3 ’ ^/2 


« . (A 2 ) 8 9. A 2 * 


4tt/j <7 3 3 


, , A(A) 3 \ ,3 , , AC A 3 * 

ii+ fTj + 4 ^(, s+ ^ - '^r 


When 8co 2 has been found, jf 2 and are determined from 

•'* ~ ft ~ l " 47T/3 a, 

r.(31). 

ya_ A!o_MJL 
/2 ( 7 g »® 4 *pCy»J 

A consideration of these formulae shows that it is immaterial what definition is 
adopted for any one of the harmonics, provided, of course, that the same definition is 
maintained throughout. 

In order to evaluate A 0 , we must eliminate 3 ©/. 


Since %* = $/<&/, 33 / 


we see that 


Hence 


®, B- = ««a U d 


dv 0 clv 0 


™ __ sin 3 /3 

dv 0 '* f,t dv 0 cos fJ cos 7 : 


sin 2 j3 \ 33/ 
cos /3 cos 7 / ' 


30/ = 33 ; 


1 (-ill _[_ ji s m iPi . x '^±(_p(y_ - 1 /<gr s S I 1M S s\2 I s ^ n2 ft §9/ (p*)l /oo\ 

4 >! +3i & + 4 <f>/ i + 2iSiPi} 5 cos /3 cos 7 8/ <f>/ - (3 


If for brevity we denote this last expression by \i, s], we have 
A 0 = ^3 [-g- (cr 2 ) 2 + 2{,j.] -joq + S [i, 6“] ... I 


*-<**+»*>^ • • • • (32) - 

ft/ = (% - a,*) & .. 

We have now the complete analytical expressions necessary for the solution of the 
problem. 


2 O 2 
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PART II. 


Numerical Calculation. 


§ 15. Determination of Certain Integrals. 
The integrals < 0 ’, pi, <£/, depend on certain others, namely— 


„ f**sin 8 "0 

Jit — .— dO 

im Jo A 8 "‘Z 




y 


(33). 


J 


After a large part of the work had been done, 1 found that these integrals tend to 
give the required results in the form of the difference between two large numbers, 
and that it would have been more advantageous to consider the integrals 


A? 


n 2n 

m 


C*" cos 2tt 0 } 

Jo 


l 


. . (34). 


It will be shown hereafter how the group (34) may easily be found from the group 
(33), and it may be mentioned that most of the results were determined in duplicate 
from both forms. 

I proceed then to consider the H, T integrals. 

Since Af = 1 — k 2 sin 2 y sin 2 6, T 2 = cos 2 y ~f k /2 sin 2 y cos 2 <f>, sin 2 /3 = w 2 sin 2 y, we 
have 


Ilf' 4 


■2m 


sin 3 j3 


/ it 2 n —2 r r 2 n —2 \ 

V L1 2'/r ~ li 2m~2J 


TZ --^ cot 2 y Tf“ 2 + ~ sin 2 y T| 4 J| 


(35). 


I now write 


J 0 

F' 


**■ cW 

A ’ 


d<j> 


p (r) = f 


0X /(1- «' 8 COS 8 (f) 

ydd 

A ’ 


E — ( A dO 
J 0 

J 'i'/r 

VO 

0 

E(y) - f Add 
Jo 


k' 2 cos 2 <j>) d<f> 


(36). 
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It will be found from Legendre’s tables that for y = 69° 49'"0, k = sin 73° 54 /- 2 


log F= '4317642, 
log F' - '2047610, 
log F(y) — '2117-987, 


\ogE= '0355145 
log E' ■= '1875655 > 
log E(y) = 9-9856045 j 


By integration by parts 


n 


2 n 

0 


2Qi -- 1) 1 +jf_ 2)l _ 2 
2 n - 1 « 3 0 


T 2 n 
0 


2(^ -1) 1 +_4i! T2 «-2 

2 n - 1 /e' 3 0 


2^-3 

(2«, - 1) «•- 


2 ^ — 4 
0 


277/ 3 

(2 n - 1) /d 3 


T 2r,— 4 
0 




Now write 


(36). 


(37)- 


6 r = -g (1 + sec 2 /3 + sec 2 y), IF = (sec 2 /3 + sec 2 y + sec 2 sec 2 y). 


The values of /3 and y are 64° 23'"712, 69° 49'"0; whence log G= "8679015, 
log H ' = 1 "4678555. Also we require hereafter log II = 1 "7182664 (see § 3). 

By differentiation 


7Asin#eos# 2n cos 3 /3 cos 3 y 2 (2n — 1) G cos 3 /3 cos 2 y 2(2w — 2)77'eos 3 /8cos 3 y 2n — ?> 

eld AF l ~ A 1 3b+3 a “ A A Aj 2 * -3 A — A^tA 

Whence, by integration, 


TT° — 

XJLo.,4.9 — 


27^ ■ 


£n° 


2n~ 2 


S' n°2n-2 + 


2n - 3 

2n 


sec 2 /3 sec 2 y n 2M _ 4 . . (38). 


On writing — 1 tan y for sin y, we find that exactly the same formula holds 
good for the T’s. 

To apply this to the determination of ng, Tq, we note that 


Also 


1I°_2 = cos 2 y F + sin 2 y E, T°_ 2 = F' — sin 2 y E' . . . (39). 

n 1 = \{F-E), To = \(F'~ E>) ..... (40). 

fC K, 

From the formula) given in Cayley’s ‘ Elliptic Integrals ’ it appears that 

^-[FEM-EFW ) 

>■ - • («>• 


n° 2 = F + 


cos fi cos 7 


r° = F f + 


sm y 


cos /3 cos y 


[F’E(y)-F'F(y) + E'F( y n 


Now ng = F, IIo is given in (40), and 14j;'"' for n = 2 , 3, 4 . . . are then given 
successively by (37). 
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Again, n 3 is given by (41), and the successive H 2 ''* are given by the general 
formula (37). 

Again (38) and (39) give 

n? — Gnl — 


n2 


5 - 2 cos-- ts cos 5 7 (° X ‘ y F + 


and by successive applications of the formula (37) we find the successive values of 

nf, n‘t. 

It is convenient also to have the series of H_ 2 , T _ 3 integrals. These are to be 
found from 

n% — njj* — sin 2 ftIi 2 0 ' l+ \ T 2 A — k' 2 sin 2 yT ^ il+2 + cos 2 yTf . . (42). 

The T integrals may apparently be derived by a similar set of formulae, but since 
at each step we divide by /c /2 , a small quantity, all accuracy is rapidly dissipated. 
Although we may safely derive one series of T integrals from a preceding one, we 
cannot so derive a succession of series, and it becomes necessary to find new formulae. 
In order to determine the T integrals, consider the group of integrals 

iw cos 2» 0 


JP n 

9 m> 


r 

J o 


rf 


d<jy e 


If we write g 


t _ cos 7 ^ an $ „ _ cos 7 


a 


we find 

0 3 + ^ 


(i + ey n 


dt 


cos ft ? cos ft 

i r 

[TO _ __j 

2m COS ft COS 2 ” 1 "" 1 7 J 0 

whence, by some easy integration#, 

TTO _ ___ 7L _ 

2 2 cos ft cos 7 9 

U\ — . y.- (sec 2 ft -f sec 2 y), 

4 cos ft cos 7 v ! ’ 

o 

U 6 = - 0 . [sec 4 ft + sec 4 y + § (sec 3 ft + sec 3 y + sec 3 ft sec 3 y) + 1 ]. 

16 cos ft cos 7 L / x / // * j 


On expanding - in powers of K f we see that 


rt&n _ TT2n I TT2n+2 i 4 ;.^ J4 TT^n+i « 

2m — 6/ 2m l 2 ^ U 2m I 2 4^ ^ 2m ’ 


When m = 0 the U integrals are easily determined. 
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The relationship between the successive U integrals is clearly 


uy _ 4 

« 2 sin" y 


Tpu-2, 
u 2»l-2 


--cot 


I now write for brevity 

x — cos /3, y — cos y, z = sin y, X — - --- -, p 
It appears that we may put 


z 

+ Z 


x + y 


2.4...2m 


1 2m -f 1 ^ t 1.3 2m "1 1.2m + 3 

2 ' 2 m ~+~2 


T 2» _ i„ h-^: 2n ~ l Jiit t n Z ..± JZ i 1-R f 2L2Lt. - .. ,/t i 

o 2^ 0/1 Om 1 ^ * O * Om i O ^ * O /I I O Q™ I A K « ’ * 


2.4 2m + 2.2 n + 4 


T 2 » _ I z a _j_ ^ 1.3...2n —JLJ , 1 2» +! > 2 . 1 • 3 2n ± 1 : 2n +3 ,, n . 

13 — " ,„. yl 2n + 0/1 , _ o i-oL i a 0 + o o„ , o a :> K + o ,i o,„ , o o„, , a C( ‘i K + • • 


2xy~ 


2(1 + 4 2.4... 2n \. 0 ~22m + 2 2 ^ 2.4 2m + 2.2m + 4 '* 


yi2n H. 77 _i_ I_ _ l- 3 • • • 2m 1 j 7 

4 “ Dm/ 2{l + 2(1+7) 2.4 . .2m 1 °° 


« , 7T 1.3...-2» —1 

;o\, t + o- ■ -v—-- wiCn 




"2 


6 16,*// * 2(1 + 4 2.4...2m \ ° 

By considering in detail the cases where n — 0 , I find 

A. o — 1, a 0 — 1 ; 

1 + 1 + phf^ 1 ~ 2p 2 f > 4 = 1 — 2 ^; 

c '° = 7 + 4 + 1(44 + .^ + ~ 3 ) + 1 


= 1 + 3 ^a - 2 p + V) + ^(1 ™ 4 /> + 7 V - W + t/> 4 ), 

C 0 = 1 %X -|- ^-X 2 . 

By some rather tedious analysis, it may be proved, by considering the manner in 
which each T is derivable from the preceding ones, that 

A _ 7 /1.3... 2m — 3 , , \ , 1.3... 2m-3 

A2n ~ 1 -2 p \2~. 4.7. 2m - 2 + + 217.72^^2 P ’ 


a 2n - l + X + r -A7_(i 


m \ 

2 m - 1 a2n ~V’ 


Ik 


1 

1 — 2 p 


P ?J I^2n-2 + 


2(i-7 2 


^ 2«-2 + 


1.3... 2m --3 

2 .4 ... 2 n — 2 


p( l — p) 


b 


2 n 



Xf a 2n 


2 m \ 2 j 

2n — 1 - 
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a. 


C 2ii 


1 


1- 2pL 
1 


9 r( , 4(1 -p)% , 1.3...2»-3 „ 7 

■ p' C 2)l _ 2 H-^ B 2n _ 2 + 2747772» - 2 p 0 ~ ^ 


1 - 2X 


(1 - X) 2 b 2n — 


3*3 

2 w\ 2 
2»-1 


C 2 


41—2 


By successive applications, starting from the values for n = 0, I find 
^0 = !» A = P) A = 9.P l 1 + P)> A — ^74 P (1 + p + Ip 2 )) 


A 

As 


1.3.5 

2.4.6 

1... 9 

2.7,10 


p (1 + p + fp 2 + fp 3 ), A 


10 — 27778 P 0 + P + 7p 3 + 7p S + ¥5p 4 )) 


p (1 + p + fp 2 + ffp S + TtP 1 + IDiP’’)' 


05-q — 1, U '2 = 1 * 4 “ X 3 cb^ —* 1 4“ X -f" i|X 2 ? — 1 4” X + 5 X 2 + I X s , 

a s = 1 + X + yX 2 + yX 3 4" ap5^ 4 )> 

a 10 — 1 + X + fX 2 + §X° 4” iftrX 4 4“ (fs^ 5 ? 


A = 1 + + 73 , A = p + A — Ip 


1 + p + ”2 (1 — Ip) 


II 


1.3 

2.4 


8 


1 + P + sP 2 + 3^5 0 + P — P 3 ) 


D 1-3.5 

^ = 2.UP 

R - 1 -' 7 n 
^ 10 “ 2 ... 8 p 


12 


1 + P + 5P 2 + fp° + K o (1 + P + Ip 2 — P”) 


1 + p + fp 2 + Ip 3 + W5"p 4 + H 2 (1 +p + 5p" Ip 8 “ fp 4 ) 


7^ 


5 


13 


1 ...9 
2~ 10 


20 


V* 4 1 


1 + P + Ip 2 + 9 P 3 + 2“1 P 4 + WsP r ' + 9^7 ( 1 + P + 7P 2 + Ip 3 — 15 P 4 — 7/ 

&0 = 1 - V, 6, = 1 + X - X 3 , & 4 = 1 + X + -|X 3 - X s , 

& 8 = 1 + X + |X 3 -|X 3 - iW b 8 = 1 4- X + -f-x 3 + |X 3 - li-x 4 - fA 

Co 


1 . 1 


^ 2/ 4 3 


+ A 4-1 (-y 3 + ^ + , 7 ) + !) 


C 2 ~ P + 3a ,3 + 3x y + 77 ) 


<7 — 1 


C 


% = t p 

__ 1_. 3 
'« 2.4 f> 


1 + P + i 0 “ 3P) + ~4J 


1 + p + Ip 2 + o (I- + p — [>') + g : (1 — ip + ip 2 ) 
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27176 P 1 '7 P + fp 2 + Ip 3 + g-^3 (l 4- p + ip 2 — p 3 ) + ! -, i( 1 +p —2p 3 +Ip 3 
2 3 P 1 + p + fp 2 + Ip 3 + 3~5P 4 + 2l77 0 + P + fp 2 — ip 3 — Ip 1 ) 


+ ,.i ( 1 + p + 0 X p 3 — 2p 3 + p 4 ) , 


1 - |A + |X7 c 2 = 1 + X - 2A 3 + fX 3 , c 4 = 1 + X + O X X 3 -2X 3 +X 4 , 

1 + X + |X a - |X 3 - |X 4 + X 8 , c 8 = 1 + X + f- fX 3 - f|X 4 - fX 5 + fX«. 


By means of these formulae I then formed a table of the II, T integrals, 
corresponding to the critical Jacobian for which y = 69°49 ; *0, k = sin 78°54'*2. 

A little consideration will show that if H 2n , n 2n , Hi n . . . are a series of II integrals, 
the A integrals as defined in (34) are as follows 

a 2 °, = nL aL = - A14, Al = A 2 nL a 2 6 , = - A 3 nL &c. 

Hence by differencing the II integrals we find the A integrals, and similarly the 
differences of the T integrals give the O integrals. 

The converse is also true, and by differencing A, O we return to II, T. 

In this way I obtain a series of values of the required integrals. It may be that 
the last decimal place is erroneous in some cases, but the results given in the following 
table are sufficiently accurate for our purpose. 


Table of Logarithms of A and O Integrals. 


n. 

log A!L 2 . 

log A 0 » 

log A 2 w . 

log A 4 n . 

log A s ». 

0 

•1064412 

•4317642 

•9049643 

1-4848513 

2*1218296 

2 

9*8518424 

9-9779679 

•1696589 

•4573383 

•8543000 

4 

9-7422028 

9-8168606 

9-9187709 

•0641097 

•2747811 

6 

9-6704204 

9-7230801 

9-7896484 

9-8771928 

9-9969099 

8 

9-6167543 

9-6573681 

9•7062560 

9-7667474 

9-8443343 

10 

9-6738269 

9•6068642 

9-6454046 

9-6909753 

9-7467821 

12 

9-5380322 

9-5658708 

9-5977435 

9-6339429 

9-6770135 


n. 

log fl" 2 - 

log 

log 

log SV-. 

log n ( ;\ 

0 

9-3901374 

•2047610 

1-0302812 

1-8667636 

2-7138342 

2 

9-0333614 

9-8993680 

•7715371 

1-6492556 

2-5319082 

4 

8-8883413 

9-7729763 

•6613027 

1-5528790 

2-4473548 

6 

8-7988101 

9-6930495 

•5898109 

1-4886133 

2-3893878 

8 

8-7344878 

9-6345770 

•5366725 

1-4398972 

2-3245748 

10 

8-6844351 

9-5884555 

•4944045 

1-4005018 

2-2861489 

12 

8-6435159 

9-5503544 

•4594781 

1-3669609 

2-2533269 


VOL. CC,—A. 2 P 











290 PROFESSOR G. II. DARWIN ON THE STABILITY OF THE PEAR-SIIAPEl) 


§ 16. The Integrals cr 3 , <x 4 , £ 4 . 


In accordance with equation (14) of the “Pear-shaped Figure” the third zonal 
harmonic is defined by 


where 


S a = sin 6 (k~ sin 3 6 — q°) (q 2 — k 2 cos 3 <]>)*,/(l — /c' 3 cos 3 </>), 
q 3 = f [1 + k 3 ~ (1 ~ 1/c 3 + k 3 )>], 9 ' 3 = 1 

The numerical values for the critical Jacobian are 


tf. 


k 2 — -9231276, (f — -5746473. 

Writing p 2 — k 3 — <f, we have 

S 3 — (p 3 — k 2 cos 3 d) v /(l — cos 3 0) ( p 3 + k' 2 sin 3 r/j) v /(k: 3 + /J 3 sin 3 c/>). 

Now let 

a = p 4 , /3 = 2p 3 «; 3 + p 4 , 

a' = p 4 K 3 , ,8' = 2p 3 /cV 3 + pV 3 , 

and we have 

(N 3 ) 3 — (a — /3 cos 3 0 + y cos 4 6 — 8 cos 6 d) (a' -f- /3' sin 3 </> y' sin 4 <f> + S' sin 6 </>). 


y — k 1 -|- 2 p 3 K 3 , 8 = K 4 , 

y = K 2 K n + 2pV\ S' = K /6 , 


The numerical values of the logarithms of the coefficients are 


Let 


log « ~ 9-0843568, log a' - 9-0496186, 

log /S = 9-8835606, log /3’ = 8-7693310, 

log y = -1748006, log y = 7-9810798, 

log 8 = 9-9305236, log 8 '' = 6-6573112. 


/(A 2ll ) = «A“ i( - /8AL + yAL - SAL 1 

j (W»i) = a TiL + /3'oL + y'ftL "1" S fi-LJ 


for n = 0, 1, 2. 


The definition of cr 3 in (8) then shows that 


era 


6 525^y*{/(A o) /(a,) - /(AJ /(n 0 ) - e[/(A„)/(n s ) - /(A,)/(!!„)]}. 


7T SUP 7 


In order to find oq and £ 4 , $ 3 must be raised to the fourth power, and we now define, 
for n = 1, 2, 3, 


/(A 2a ) = ot 2 AL — 2a/3A| t + (2ay + /3 2 ) A 2 4 (l — (2a8 + 2/3y) A| t 

+ (2/38 + y“) Ajj,, — 2ySA 2 ° + S 2 A 2 f ( , 
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/(-0 2 ,) = *' 2 nl + 2a' / S'0| l + (2«y + ,Q' 2 ) O 4 , + (2a' 8' + 2,8'/) 

+ (2/3' 8' + / 2 ) .fil„ + 2/ S'< + S' 2 ft£. 

From the definitions of cr 4 , £ 4 in (8), we see that 

U = l -iJf'-- 1 { - [/(A., )/(!»«) -/(A,)/(nJ] + G [/(A s )/(li tl ) - /(A t )/(n s )] 

- //[/(A,)/(o 4 ) -/(A,)/(a)]), 


6 cos 3 /3 cos 3 y sin /3 


~ 


sin 3 7 


-/{n.)/( a,)] 

-se[/(A,)/(n 1 )-/(A t )/(n > )]}. 


The computations (which were in this case actually made from the corresponding 
formulae involving the II, T integrals) gave 

og = -0136866, £ 4 — -00009246, cr 4 = -00176135. 

These have to be used in a formula which also involves Now % ?> denotes 

or what should be the same thing, P/QJ. The formulae in the “ Pear-shaped 
Figure” with y — 69° 49 ,- 0, k — sin 73° 54'*2, give 

P^Qji = -351697, = -351744. 


Thus the two functions, which should be identical in value, differ by "000047. I 
think that if I had taken y = 69° 48'"997, k = sin 73° 54'"225 (the actual numerical 
solution for the critical Jacobian, although not fully stated in the “ Pear-shaped 
Figure ”) this small discrepancy would have been removed. However, the difference 
is quite unimportant, and as <2I 3 generally means P 1 1 Q 1 1 , I take the former value and 
put log iHg = 9-54617. 

With this value I find the required result, namely 

[I (°h) 2 + 2t,] - ior 4 = - -00050012 ..... (43). 


§ 17. The Integrals &>/, pi, <f> ; \ 

Any harmonic S*, where i and s are both even, whether in the approximate form 
of “ Harmonics ” or in the rigorous form, may be written 

Si = (a — h cos 2 6 -j- c cos 4 6 — d cos 6 6 ~f- ...)(«' + h’ sin 2 <f> + o' sin 4 <f> + d' sin 6 </> + ...). 

Each series is, of course, terminable, the number of terms in each of the two factors 
being \i -j- 1. 

For the determination of the at, p integral^ this must be multiplied by ($ 3 ) 2 . It 

2 p 2 
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can be seen, without actually writing down the product, how the coefficients will 
occur; I write therefore those coefficients as follows :— 

l 0 “ cta , 7j < 1 let ba, ™ cty — bft cot, b ™ o .8 j by eft -1- c/ot, Ac. 

m 0 — afa, rn 2 = a'ft' -j- 6V, m 4 — ay' + l/ft' + da!, m R = a! S' + b'y' -f- c'ft' + d'a', &c. 


Next let 


f {ft-lft) - 4Ao n 4-^211 H" 4^-2« 4-^-L. 


/ (o 2ll ) — w 0 n| B + m 2 nl n + w 4 o| b + w 6 o| a ... 
Then it follows from the definitions of toft p/ in (8) that 

< = {/(A„)/(%) -/(Aj/(Ho) - ff[/(A 

„._i.»A^«’p /(Ao)/(a2) - /(Aj)/(f!o)! . 


for ~ 0, 1,2. 


/(A 2 )/(fto)]} 5 


*?r siir 7 Bin p 


It is, of course, necessary to reduce the two factors of $* to the required forms. 
The harmonics of the second order are 


$/ = (k 2 sin 2 6 — p) ((/ 


' 2 -- 4“ cos 2 </j), (.v = 0, 2), 


and I find qft 2 = ‘3197540, qft = ‘9623311 ; whence we may find a, b, a!, b' for 
these harmonics. 

For the harmonics of the fourth and sixth orders I take the formulae of “ Harmonics,” 
and attributing to the parameter ft its value ‘0399726 (or more shortly ‘04 in the 
terms of the sixth order), 1 reduce (/a), C/ (<j>) to the required forms and determine 
a, b, c, &c., a', V, cj, &c. The numerical values of these coefficients are given in the 
tables of § 20 hereafter. 

It may be well to remark that p 0 is needed (but not toft), and in this case S Q = 1, so 
that a = a' == 1. 

It seems useless to go in detail through the tedious operations involved in carrying 
out this process in the several cases. 

Approximate formulae are given for the <f> ; s integrals in § 22 of “ Harmonics.” The 

f p 0 clcr of that paper is the same as -§7r/c 3 f [ <b d0 dtf> of the present one, and 

J siir p JJ 

the factor there written M is P - . Hence it follows that 

sir j3 


p/C/) 3 p dcr of “ Harmonics. 


In order to apply this to the harmonics of the second degree, it must be borne in 
mind that a different definition of Sp (s = 0, 2) is being used here. If [</>J, [</>/] be 
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the values which would be found from “ Harmonics ” without this correction, and if 
(f> ? , are the required values, it appears that 


: 




where «, e, a, e are the coefficients specified in § 12 of “ Harmonics.” 

The approximate values found in this way for all the <f> integrals are very near to 
the more correct values, and might have been adopted throughout without material 
error. But as there was not much certainty that the approximation was a good one 
—and indeed for S a was probably bad—I also found all these integrals, excepting 
<j j i 6 3 , <^ 6 4 , by the method now to be described. 

From (6) and (8) it appears that 




If, therefore, we write 


7 T siir 7 J \A,' 



CS /) 3 

AF 


dO d <[>. 


jf (A 2 ,) —- — 2abA$ n -f- (2 ac + Ir) A% n — (2 ad + 2 be) A 2h + • • • 

f{Fl 2) i) ^ (ij "HL 2 a b Fl 2n (2a c T~ b w ) (2a d -j- 2 be') Fl\ n -J- . e . 


we have 


H z:z 


t r siir 7 


[/ a) / (^o) ,/ (A 0 ) J (fl-a)]* 


The following table gives the results for all the a >/, pf, <£/ integrals :— 


Table of Logarithms of <f>, w, p Integrals. 


t s. 

log 

Approximate log c/>/ 
from “ Harmonics.” 

log (Uj*+ 10. 

log pi s + 10 . 

0 

0 

— 

— 

— 

7-63099 

2 

0 

9-00516-J0 

9-00515-10 

7-67371 

7-02716 

2 

2 

7-03970-10 

7-03973- 10 

(-) 5-68193 

(-) 5-05256 

4 

0 

9*69080- 10 

9*69323 - 10 

8-03358 

7-35625 

4 

2 

1-72664 

1-72729 

(-) 8-33367 

(-) 7-59132 

4 

4 

3-81541 

3-81612 

8-29058 

7•37446 

6 

0 

9-71219-10 

9-69305- 10 

7■97301 

7-32602 

6 

2 

— 

2-20562 

(-) 8-72778 

(-) 7•94094 

6 

4 


5-29999 

9•10091 

8-13161 
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§18. The Integrals B/. 


Adopting the notation of the last section we have 

f / (v) = a + h (v 3 - 1) + o (jd - l) 3 + d (y 2 - I) 3 + . . . 

Let v — —. - , and k sin il> — sin v, so that at the surface of the ellipsoid where 

K Sill ^ ' A x 

ij/ = y, we have y = /3. 

Then 

P/ (c) = a -f- h cot 2 y “h c cot 4 y + . . . 

P/ (*» 0 ) = a + & cot 2 /f f c cot 4 /3 + . . . 

P” civ 


Now 


and 




1-^ 


city 


-s 1 + ' K \ ^/ Ji _tyz surety) 


k \ sec x d\jj 


civ 


I* - Li/ Ml \ (§ff{ty> - 


,2 _ L+J?Y 

1 - /3/ 


Hence 


gj/ — k (a + h cot 3 /3 + c cot 4 /3 + . . . ) 3 | 


sec % city 


o [a -j~ & cot 3 % + c cot 4 . .] 2 


We have, in § 4 of the “ Pear-shaped Figure,” the rigorous expression of this 
integral for harmonics of the second order, viz. :— 


k (1 — 2q?) (1 — q? sin 3 y) 3 [ F( y) E (y) sin y cos y cos fi 
2 q?q', % sin 4 y [ q/ q/q'f q' s $ (1 — sin 3 y) 


s = 0, 2. 


The values of < 2 o 3 , g 1 / have been already given, and thus all the quantities involved 
are known. 

The two factors of (viz., fp and are given in approximate forms in 

“ Harmonics,” and therefore, if we made allowance for the different definition of © ; * 
adopted in that paper, we might calculate The .computations I made showed 

that the results obtained in that way would have been sufficiently exact, but as it 
was clear that the approximation to the functions was not very close, and as the 
computation is tedious, it seemed better to find the by quadratures. 

In order to do this I divided y or 69° 49' by 12, and took 5° 49'-^- as the common 
difference, say 8. I then computed sec y, a + b cot 2 y + c cot 4 X + • • ■ > and 
sec y 4- (a + h cot 2 y + c cot 4 y . . .) 2 for values of \jj — 0, 8, 28 , . . . 128 or y. 

As a fact the first five or six values need not be computed because the early values 
of the functions to be integrated are practically zero. The ordinary formula; of 
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quadratures are inappropriate for these integrations, because the function, say u H , to 
be integrated increases so very rapidly. I therefore take an empirical and integrable 
function, say v n , which is such that v u = zq 2 , v n = u u ; the quadratures may then be 
applied to u n — v m and the result applied as a correction to | v d\p, In fact this 
correction is always very small, and we might well be satisfied to use J v dijj , which is 
very easy to calculate. 

The empirical function v is given by 


(*-7) loggia 

v = u u e 8 u n . 


Then when \fj = y, v = u n ; when \fj = y — S, v = u u ; and 



In all the cases I have to consider the exponential term is negligible, and the 

^ 12 ^. 

integral is lo ^ . 

For the quadratures we have 


V 12 — Uj2, 5 'V 


10 


u 


13 




> v 9 = u i*[—) > &c - 


u. 


U2. 


and the equidistant values of the function, to be integrated (arranged backwards), are 


u- 


12 , 


U- 


12 , 


A A / 7 fil \ 2 MlV MF 4 

0, 0, U 10 ~ U l2 ( — ) , U< 3 — U 12 ( ) , U s — U 12 1 


"12 


, &c. 


The first two are zero, the next three or four are found to be sensible, and the rest 
are insensible ; hence the quadrature is very easy. 

The 3 Si integrals are found thus :— 

_ gr. __ 1 . d an . /_, \ _ 2 ®i‘ h + 2c cot 2 # + 3d cot* £ + ... 

; 1 (v 0 ) dv 0 Vi Vo/ sin p ci + b cot 3 p + c cot 4 p + ... ' 


The following table gives the 2C, integrals • 
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Table of Logarithms of the ( Jf., JJ integrals. 


i. 


log fi/ + 10 , 

log HA 

2 

0 

9*69312 

*09295 

2 

2 

9 * 33300 

•40665 

3 

0 

9*54617 

•20467 

4 

0 

9*44928 

*28206 

4 

2 

9*25987 

*41239 

4 

4 • 

9*06489 

: 

-44858 

6 

0 

9*24383 

| 

-35876 

8 

2 

9-16199 

•41249 

6 

4 

9-02369 

-44195 


§ 19. Synthesis of Numerical Results; Stability of the Pear. 

In the following tables and remarks I collect together some of the results which 
occur in the course of the work. The final places of decimals as given have, perhaps, 
in many cases but little significance :— 


i. s. 

(i.) 

a 8 - fb. 

(2.) 

log ($3 - % s ) <l>i S = log Of. 

(3.) 

(4.) 

/,» sm 2 fi \ 

1 3 2 co s/3 cosy/^ * 

(3) + (4). 

Bf. 

2 0 

2 2 

- -141617 
•136417 

(-) 8-1562926- 10 
6-1745705 - 10 

-0029865 
- *0000247 

- -0011976 
•0000127 

•0017889 
- •0000120 

4 0 

4 2 

4 4 

•07033 

•16978 

•23558 

8-53794 - 10 

•95653 

3-18755 

*005214 
- -008965 
*005594 

- -002555 
•004390 

- -002664 

*002659 
- -004575 
•002930 

6 0 

6 2 

6 4 

•17638 

*20649 

•24609 

8-95864 - 10 

1-52052 

4-69108 

*004067 
- *019041 
•032054 

- *002384 
•009820 

- -015232 

•001684 
- -009221 
•016822 


For all harmonics higher than those of the second degree — ft/ is the coefficient 
of stability. Since in all these cases this expression is positive, the ellipsoid is stable 
for all such deformations. 

If U -j- 8 U be the energy function for the pear, whose variations for constant 
moment of momentum are considered by M. Poincare, we have in our notation 

U+'bUA= - 11 + | (Aj - A,) (<P + Swf 
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It is easy to show from our analysis that for the deformation f%$ 2 , 

or C" 


O TLf 2 

8£r= W (/i) M <«»-«»)* + 


2irp fl 


and that the corresponding expression with ft 3 in place of r 3 holds good for the 
deformation f 2 S. z . 

Forestalling the results obtained below, it may be stated that for f z S. z 


and for f z S z 


{-•01433 + -03959} ; 

^ V Q 

W = ~ (fif {-00015 + -00002}. 


Thus in both cases 3Z7 is positive, and this shows that the Jacobian ellipsoid is also 
stable for the ellipsoidal deformations. The fact, that 8 E (the variation of my 
function of energy for constant angular velocity) is negative for the deformation S. 2> 
illustrates the truth of M. Poincare's remark (‘Acta Math.,’ 7, p. 365): “Si au 
contraire la rotation de la masse fluide etait determinde par celle d’un axe rigide 
(comme dans les experiences de Plateau par exemple), tout deplacement produirait 
une resistance passive et 1’ellipso'ide de Jacobi serait toujours instable.” 

I have in (32) written 



+ iWpif 


sin 3 /3 


4 cos /3 cos 7 



The following table then gives further stages in the work : —■ 


i. 

s. 

[*» *•] 


Bf/Ci*. 

2 

0 

•00014032 

- *00022329 

- *12482 

2 

2 

•00000072 

+-00000097 

- *08059 

4 

0 

•00009937 

•00020486 

•07705 

4 

2 

•00000276 

•00000231 

- *000506 

4 

4 

•00000001 

•00000001 

•0000019 

6 

0 

•00002835 

•00003118 

•01852 

6 

2 

*00000403 

•00000256 

- -000278 

6 

4 

•00000002 

•00000001 

• 00000034 




•00024190 -*00022329 




•00027558 

+ *00024190 


% [!K> 2 + 2fd-^ 4 = 

- *00050012 

( P -A2 




-- 

= + -00001861 



Aq — 

- *00022454 




i I! 

+*00001861 



Numerator. 

- *00020593 




vol. cc.— a. 2 Q 
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The next step is to find 
r sin /3 


cos 3 /3 


eos/3 cos 7 \ D + /ey ’ 


sin ft / cos 3 /3 '■ 
4 Dy 3 2 cos /j cos y \ D — /cy 


where Z) 2 = 1 — kV 2 . 

The numerical values are log D = 9‘9840165, log L = ‘6454565, log M = ‘9591960. 
From’these we obtain c = L<f> 2 , tl = M(j) 2 ; whence 


A 

<h c ~ 

- -055837 

A 2 

•000804 

I) = 

•031701 

Denominator = 

- -023332 


In accordance with (32) the Numerator divided by the Denominator is 


and I thus find 


I °Si^ = <-)!''94578. 


It was found in § 7 of the “ Pear-shaped Figure ” that the angular velocity of the 
critical Jacobian was given by -—= T4200. Accordingly, the square of the angular 
velocity of the pear being w 3 + 8 « 2 , we have 

w 2 + 8<u 2 = ft) 3 [l - -124314c 2 ]. 

From the formulse (31) I then find 

— ‘15068e 3 , / 3 2 = ‘50839c 2 . 

The other f 3 are equal to ~ e 3 , and are given in the preceding table. From (28) 
and the definitions of a, f>, C, ft it appears that the moment of inertia of the pear is 


Aj — A r 


2? rpJc 0 [_ ^ a 'a-* 2 a ^ 3 


With log a = 9‘8559758, I find 

^-^=IS [1+ ' 131011e2] - 

The angular velocity of the pear is 

y/ (ft / 2 + 8 ft) 2 ) = w [1 — -062157e 2 ]. 
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Multiplying these last two expressions together, we have the moment of 
momentum of the pear ; it is 


2 TTpJc 0 


[1 + •068854c 3 ]. 


It follows that, whilst the angular velocity of the pear is less than that of the 
critical Jacobian, the moment of momentum is greater. This result would afford a 
rigorous proof of the stability of the pear if the numbers were based on a complete 
solution of the problem. But as we have not determined an infinite series of new 
harmonic terms, it becomes necessary to consider how the result might differ if the 
hitherto uncomputed terms were added. 

tc . j - a .— j-'L„ ...__ j-i — ,.j? u. • _•_v frv i (-®/) 2 l -1 ;.c 


If e denotes the uncomputed portion of the infinite series % j \i, s] -j—j, and if 

A denotes the addition to be made on that account to any of the results as already 
computed, we have 


A \4ttp) ~ ‘023332 ’ &nd A \rf t 


•023332 x 142 * 


Whence 


Since 


A [V(co 3 + Sco 2 )] = |coA “ = co [301*8346ee 2 ]. 


/*_^4 2 | 

2 ~ ^ 4 t rp ( 7 3 


f-2 _ ?i_ .,2 _ M JL 

H C* 4 nrp C* 


A ^ = 4 a S = - 10 '” 

A 2 _ _ JL a .^ Ct)2 _ _ 1 a3’4563848 

A/2 - C* A 47T P “ 1U 


Then 


A G, - A .)-= %£ [‘ A/, - | A//] = [- 833-7892 + 39-7472] ce* 


3Jf 2 a 


(- 794-0420) ee 3 . 


Therefore 


^(w 3 + Sw 3 ) = [1 — -0621568e 3 + 301'8346ee 3 ] 

A- - a 4, = —~-[1 + -13101068c 3 - 794-0420ee 3 ] 

Airpk 0 

By multiplication we find that the moment of momentum is 


3 M*aco 

2irpk 0 


[1 + -0688539c 3 — 492‘2074e 3 e]. 
2 Q 2 






300 PROFESSOR G. H. DARWIN ON THE STABILITY OF THE PEAR-SHAPED 


The coefficient of e 2 is positive and the pear is stable, provided that 

492'2074e < ‘0688539, 
or e <C '00014. 


Inspection of the table of numerical results shows that the zonal harmonic terms 
contribute by very far the larger portion of the sum. Now the sixth zonal term was 

[6, 0] + = '00002835 + '00003118 = '00005953. 


This is about f of the critical total '00014. The pear is then stable unless the 
residue of the apparently highly convergent series shall amount to 2-j times the 
contribution of the sixth zonal term. Such an hypothesis appears profoundly 
improbable, but I have thought it expedient to make a rough determination of the 
contribution of the eighth zonal harmonic to the sum. 

If we take k as equal to unity, S s — (/x) Cg (< f> ) = P 8 (,u), and we easily see that 

the formulae (8) reduce to 


Ct) s 


cos 3 y 


ir /•■Iff (1 q. gi n 2 y gi n 3 CO g Q 


Jo Jo 


P 8 — 


(1 — sin 2 7 sin 3 6 J 
6 cos 2 7 f ^ cos 6 


(S^S 8 d0dcl>, 


w 


ir siny Jo Jo 1 — sin 3 7sin 3 6 


(s s ys a ddd<f>. 


In these integrals (fi only enters through (S 3 ) 2 or [|J 3 (yx)] 2 [C 3 (<£)] 2 . 
Now 

T[C = 

J 0 


1 Q 1 U K 

«' + F + ^y + ^8' 


Hence 


w 8 — I K cos 2 y | 
3/Ycos 3 7 r** 


2.4 ; 1 2.4.6 

\trK, where K — *1452. 


[fa (/»)]*-P 8 (fi) d,o 

0 (1 — sin 3 7 sm 3 6 y v 0 u 1 


cos 6 


Ta Cfa 


sin 7 J o 1 — sin 3 7 sin 3 0 

In these integrals 

Pg (fx) = T |s [6435 sin 8 d - 12012 sin 0 d + 6930 sin 4 6 - 1260 sin 2 d + 35] 


DPs (m)] 3 = « — y3 cos 2 d + y cos 4 d — 8 cos 0 d, 

where a, yQ, y, 8 have known numerical values. 

The integrations may of course be effected rigorously, but it seemed far easier to 
determine them by quadratures. I therefore computed the values of the functions to 
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be integrated for-f? — 0, 15°, 30°, 45°, 60°, 75°, 90°, drew curves on squared paper, 
and counted the squares on the positive and negative sides of the axis. 

In this way I find logwg = 6*765, log p s = 5-653. 

The integral </> g is found at once from § 22 of “ Harmonics ” with /3 — 0. This 
gives cf> s = tV, or log </> 8 = 9-247. 

If P 8 (ju,) he expressed in terms of cosines of 0 we have 

P 3 (/x) = a — b cos 3 6 + c cos 4 0 — cl cos 6 0 + e cos 8 0, 

where a = 1, b = 18, c = 74-25, d — 107-25, e — 50-273. 

Then we may, as in §18, put 

u„ — P s (v) = a + b cot 2 x + c cot' 1 x + d cot 0 x + e cot 8 y- 

As was done in that section, I then computed u n and u n , and so found the integral 
of the empirical function. The result gave 

log & 8 = 9*191 ; whence log3B 8 = "370. 

It may be admitted that the determination of ^ 8 , 1$ 8 is not wholly consistent with 
that of the previous integrals, since I only assume k to be unity in as far as the values 
of a, b, c, d, e are affected. 

Applying these values as before, I find & 3 —& 3 = T97, log O 8 =8"540, B s = -000092, 

yr = "0027, and 
C'8 

[8, 0] == *00000051, ® = -00000025. 

Hence that part of e (the uncomputed residue of the series) which depends on the 
eight zonal harmonic is only about "0000008. The contribution is so insignificant 
compared with the critical total "00014, that I have not thought it worth while to 
make estimates for the tenth and twelfth harmonics. 

It may then be confidently asserted that the pear is stable. 

In the course of this estimate we have also found f s = — '0027e~. 


§ 20. Second Approximation to the Form of the Pear. 

Extracting the numerical values ot the f’s from our results, we find that the 
inequality of the critical Jacobian ellipsoid is 

eS. d + e■■ [-15068 S z + '50839 Sf + *07705 — *000506 + *00000019 Sf 

+ -01852 S 6 - -000278 Sf + ‘00000034 £ 0 4 - ? S a 6 + '0027 S s — . . .]. 
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In order to give this expression a clear meaning, it is well to define the several S’ s. 
S s = (k 9 sin 3 0 — <f sin 0) (<j' 9 — k' 9 cos 3 <f) ^/(l — k' 9 cos 3 <f>), 
where k 2 = -923128, q 9 = '574647 

/c' 3 = -076872, q /3 = -425353. 

For the other harmonics we have 

S* = (a — b cos 3 6 + c cos 4 6 — e7cos 6 6 + . . .) (a'4- V sin 3 sin 4 sin 6 . . .), 

where the values of a, b, &c., a\ V, &c. are as given in the following table :—• 


i. 

s. 

a. 

b. 

c. 

d. 

e. 

2 

0 

•603374 

•923128 




2 

2 

- -039203 

•923128 




4 

0 

1 

5-450 

4-901 



4 

2 

- 1-7988 

36-006 

44-805 



4 

4 

•0839 

-7-975 

95-574 



6 

0 

1 

12*6 

29-984 

18-834 


6 

2 

-8-4 

121-8 

439-425 

320-523 


6 

4 

3-78 

-338-312 

; 3680-303 

4482-844 


8 

0 

1 

18 

74-25 

107-25 

50-273 


i. 

s. 

a'. 

V. 

c\ 

d'. 

e. 

2 

0 

•603374 

•076872 




2 

2 

- *039203 

•076872 




4 

0 

•8036 

•3718 

•0280 



4 

2 

-1*0666 

1-8135 

•1865 



4 

4 

1-0136 

-8-0266 

8* 



6 

0 

•5989 

-6888 

•1512 

•0000 


6 

2 

-1-1408 

1-5349 

•4404 

•0264 


6 

4 

1-0305 

-7*704 

6-944 

•704 


. 8 

0 

1 

0 

0 

0 

0 


The surface of the pear is determined by measuring a certain length along the arc 
of curves orthogonal to the surface of the ellipsoid. By equation (22) it appears 
that that length measured in the direction of the positive normal is 
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f>» [«&, + tfiS- + $e? (i + ~ - no) (S,y 

In order to construct a figure it will be convenient to adopt as unit ol length c, 
the greatest axis of the ellipsoid which is deformed. We know that 

c = - , b = k cot B, a = ^ cos ^ , so that b = c cos B, a = c cos y, and the mass or 

sin ft sm ft ' 

3 cos /3 cos 7 


the ellipsoid is f 7 rpc 3 cos /3 cos y. But since the mass of the pear is inpk 
where k 0 s — F J (1 + <? 3 cr 2 ), it follows that it is 


sin 3 /3 


%Trpc s cos /3 cos y (1 + - 0136866e 3 ). 


Hence the mass of the pear is a little greater than that of the ellipsoid whose 
deformations we shall draw, and the protuberances above the surface slightly exceed 
in volume the depressions below it. 


We have p 0 


c cos B cos 7 


c cos B cos 7 




(i 


sin 2 ft sin 2 6 f (cos 2 7 + K n sin 2 7 cos 2 (f>y 


and the expression for the orthogonal arc, measured from the ellipsoid to the pear, is 
therefore 


Po 


eS 3 + 


c J ( 2(1 — sin 2 B sin 2 6 ) 


+ 


2 (cos 2 7 + k ' 2 sin 2 6 cos 2 (f>) 

— |(1 + sec 3 B + sec 2 y) \ + ZfiS? 


It appears to me that it will afford a sufficient idea of the corrected form of surface 
if I draw two principal sections, namely, first, a section through the axis of rotation 
and the longest axis of the ellipsoid, and, secondly, a section at right angles to the 
ay is of rotation. It is not worth while to consider the third section drawn through 
the axis of rotation and the mean axis of the ellipsoid, since it will hardly differ 
sensibly from the uppermost figure shown in the “ Pear-shaped Figure.” 

For the sake of brevity I will call the first and second sections the meridian and 
the equator. 

The three ellipsoidal co-ordinates v, 6, (j) of any point are connected with x, y, z by 
the relationships 

x — c sin y. (/cV- — 1)* (1 — /c 3 sin 3 0)* cos cf>, 
y = c sin y . k (v 2 — 1 )* cos 9 sin </>, 
z — c sin y. kv sin 0 ( 1 — k % cos 3 <p)K 

The equation to the surface of the ellipsoid is v = —~— = . 1 ;, • 

1 1 k sm 7 sm» 
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The equation to the meridian plane in rectangular co-ordinates is simply y — 0, 
that to the equator is x = 0. 

In ellipsoidal co-ordinates the equation to the equator is simply <j> — \tt, hut the 
equation to the meridian is peculiar, for it is in part represented by 0 = and in 
part by <j> = 0. 

The curve 0 = <f> = 0, which defines the limit between the two regions where 

the equation to the plane has different forms, is clearly the hyperbola 


In the region from z — oo and * small down to this hyperbola the equation is 
0 = y7T; and between the origin and the hyperbola it is <{> = 0. 

If we follow the arc of the ellipse from the extremity of the c axis we begin with 
0 — T?r, (f) = 1-77, and 0 remains constant whilst <f> falls to zero. Then cj> maintains a 
constant zero value whilst 0 falls from hr to zero. 

On the side of the origin where z is negative, 0 is of course negative and undergoes 
parallel changes. 

The hyperbola 0 = \ 2 tt, cf> = 0 cuts the ellipsoid so near to the extremities of the c 
axis that an adequate idea of the deformation may be derived from the two extreme 
values of <f>, namely, \tt and 0. I have also thought it sufficient to compute the 
deformations for 0 = 0, 30°, 60°, 90°. We thus obtain the following scheme of values 
of 0, <jf>, together with the corresponding rectangular co-ordinates (with c taken as 
unity), at which to compute the deformation :— 


Meridian (y — 0). 


0 = 

0 

o 

o 

<t> = 

90°; 

; Z = 

1, 

X = 

0 

e- 

so 

o 

o 


0; 

z = 

•961, 

X = 

•096 

e~ 

o 

o 

<£ = 

0; 

z = 

•832, 

X = 

T91 

0 = 

30°, 


0; 

Z — 

oo 

o 

X = 

•303 

0 = 

0, 


0; 

z = 

0, 

X = 

•345 


Equator (« = 0). 


e- 

- 90°, 

cj) = 

= 90°; 

z=l, 

y — 0 

0 = 

r 60°, 

$ = 

= 90°; 

z = *866, 

y= -216 

0 = 

= 30°, 

(j> = 

O 

O 

n 

z = *5, 

y = ‘374 

e- 

= o°, 

4> = 

= 90°; 

z — 0, 

y- -432 


It did not seem to be worth while to compute the deformations due to the eighth 
zonal harmonic, since it would be quite impossible to show them on a drawing of any 
reasonable scale. 

In order to exhibit the magnitudes of the contributions of the harmonics of the 
several orders, I give the normal departures 8n at the points z = ± 1, x — 0, y = 0. 
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Term of first order . . . 

S s 

± •148227c 


Terms of second order pro- 




portional to e 2 . . . . 

W 


- -010986 


8 a 

•061844 



Sf 


•000751 


s 4 

•092715 



SJ 

•000849 



S4 4 

•000000 



s 6 

•026647 



s<? 

•002011 



s 6 4 

•000001 




• 184067 

- -011737 



- -011737 




• 172330e 2 



The following are then the results for the normal departures at the several points 


whose rectangular co-ordinates are specified :— 


Meridian (y — 0). 


i h * 

—h '961, x 
—h '832, x 
± '480, x 
0, x 


0, Sn= ± T482e + '1723c 3 . 
•096, Sn = ± '0932e + '0858c 2 . 
•191, 8 n - ± '0189c + '0103c 2 . 
•303, Bn = '0223c - '0033e 2 


*5}/I ^ 


i * s\ r\ a r* 9 . 


Equator (^ = 0). 

z = ± 1, y = o, Sft = ± '1482c + -1723c 2 

z = ± -866, y = -216, Sn = ± -0300c + -1265e 3 . 

z = ± -5, y = -374, 8n = qp -0354c - ’0220c 2 . 

* = 0, y = -432, Sft = - -0095c 2 . 

In order to draw a figure I take e — |. Throughout most of the arc of the 
ellipsoid the approximation is probably good, but at the vertices, which are just the 
points of most interest, it is pretty clear that we are using a somewhat extreme value 
for e. The results are • 

2 n 


von. cc. —A. 
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Meridian (?/ = 

0 ). 




Equator (x — 0). 



z = 

i, 

x — 0, 

Sn = 

+ 

•117. 

z = 

i, y = o, 

Sri = 

+ T17. 

z = 

•96, 

x = -096, 

S n — 

+ 

•068. 

0 = 

r-H 

qr 

ii 

OO 

Sw — 

+ -047. 

z = 

•83, 

x — T9, 

Sn = 

+ 

•012. 

z = 

•5, y = -374, 

Sn = 

- ‘014. 

z ~ 

•48, 

8 

II 

CO 

o 

Sn = 

— 

•on. 

z = 

0, y — -432, 

Sn = 

— -002. 

z = 

o, 

* = -345, 

Sn = 

+ 

•001. 

z = — 

•5, y = -374, 

Sn = 

+ -003. 

z = 

- '48, 

x = ‘30, 

Sri = 

+ 

•010. 

z = — 

’866, y = ‘216, 

Sn = 

+ -017. 

z == 

- '83, 

x — T9, 

Sn = 

— 

•007. 

z = — 

1, 27 — 0, 

Sn = 

— -031. 

% rzz 

- -96, 

x = •096, 

Sn = 

— 

•025. 

N.ll- 

-For z = fit: '866, 

Sn is 

in both 

z ~~ 

- 1 , 

x —■ 0, 

Sn — 

— 

•031. 


cases positive . 



These numbers are set out graphically in the annexed figure. It will be noticed 
that whereas the protuberance at the positive end of the z axis is great, the 


B 



A 



deficiency at the negative end is almost filled up. We may describe the general 
effect by saying that the Jacobian ellipsoid is very little changed, excepting at one 
end of its longest axis, where it shoots forth a protuberance. 


Summary. 

If a mass of liquid be rotating like a rigid body with uniform angular velocity, the 
determination of the figure of equilibrium may be treated as a statical problem, if 
the mass be subjected to a rotation potential. 

The energy, say W, lost in the concentration of a body from a condition of infinite 
dispersion is equal to the potential of the body in its final configuration at the 
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position of each molecule, multiplied by the mass of the molecule and summed 
throughout the body. In the system, as rendered statical, it is necessary to add the 
rotation-potential to the gravitation potential before effecting the summation. That 
portion, say T, of the whole lost energy which arises from the rotation-potential is 
simply the same thing as the kinetic energy of the mass, when the system is regarded 
as a dynamical one. If we replace W + T by E to denote the whole lost energy of 
the statical system, the condition that the surface shall be in equilibrium is that the 
variations of E for constant angular velocity shall be stationary. E must then be a 
maximum or a minimum, or a maximum for some variations and a minimum for 
others. 

It might appear at first sight that the condition for the secular stability of the 
figure is that E should be a maximum for all variations, and this is so if certain 
constraints are introduced ; but in the absence of such constraints the figure may be 
stable although E is a minimax. 

It has been shown by M. PoincakIs that the stability must be determined from the 
variations, subject to constancy of angular momentum, of the total energy of the 
system, both kinetic and potential. The two portions of the total energy, say Z7, are 
again W and T ; but whereas E involves the lost energy W of the system under the 
action of the gravitation potential, U involves the potential energy which is equal to 
- IF. Thus U is equal to — W + T. 

The variation of U with constant angular momentum leads to results for the 
determination of the figure identical with those found from the variation of E with 
constant angular velocity. But there is this important difference, that to insure 
secular stability U must be an absolute minimum. It appears, in fact, that, in the 
case of the pear-shaped figure, while E is actually a maximum for all the deforma¬ 
tions but one, it is a minimum for that one, which consists of an ellipsoidal strain of 
the critical Jacobian ellipsoid from which the pear-shaped figures bifurcate (§ 19). 

But M. Poixcare has adduced another consideration which enables us to determine 
the stability of the pear by means of the function E, without a direct proof that U is 
a minimum for all variations. For he has shown that if for given angular momentum 
slightly less than that of the critical Jacobian ellipsoid, the only possible - figure 
is the Jacobian, and if for slightly greater angular momentum there are two figures 
(namely, the Jacobian and the pear # ), then exchange of stability between the two 
series must occur at the bifurcation. If, on the other hand, the smaller momentum 
corresponds with the two figures and the larger with only one, one of the two 
coalescent series must be stable and the other unstable. Now it has been proved 
that the less elongated Jacobian ellipsoids are stable, so that if the first alternative 
holds the stability must pass from the Jacobian series to the pear series; and if the 
second alternative holds the pear series must be unstable throughout. The question 

* For the sake of simplicity we may speak of a single pear, instead of two similar pears in azimuths 
180° apart. 


2 E 2 
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of stability is then completely determined by means of the angular momentum of the 
pear ; if it is greater than that of the critical Jacobian the pear is stable, and, 
if less, unstable. 

It suffices then to determine the figure by means of the variations of E with 
constant angular velocity, and afterwards to evaluate the angular momentum. 

It was proved by M. Poincar^, and repeated by me in my previous paper, that 
the first approximation to the pear-shaped figure is given by the third zonal 
harmonic inequality of the critical Jacobian ellipsoid—zonal with respect to its longest 
axis. In proceeding to the higher approximation I suppose that the amplitude of 
the third zonal harmonic is measured by a parameter e, which is to be regarded as a 
quantity of the first order. We must now also suppose the ellipsoid to be deformed 
by all and any other harmonics, but with amplitudes of order e 2 In the first 
approximation the lost energy IF is proportional to e 2 , but it now becomes necessary 
to determine W as far as the order e 4 . A change in the sign of e means that the 
figure of equilibrium is rotated in azimuth through 180°. Such a rotation cannot 
affect the value of the energy, and it thus becomes obvious that the odd powers of e 
must be absent from the expression for IF. We have further to find the moment of 
inertia of the body as far as the terms of order c 2 , and thence to find the kinetic 
energy 1\ The function E is equal to IF + T. 

In order to attain the requisite degree of accuracy, it is convenient to regard the 
pear as being built up in an artificial manner. I construct an ellipsoid similar to and 
concentric with the critical Jacobian, and therefore itself possessing the same 
character. The size of this new Jacobian, which I call J, is undefined, and is subject 
only to the condition that it shall be large enough to enclose the whole pear. The 
regions between J and the pear being called R, I suppose the pear to consist of 
positive density throughout J and negative density throughout R (§ 1). 

The lost energy of the pear consists of that of J with itself, say \JJ ; of J with 
R , which is filled with negative density, say — JR; and of R with itself, say | /tA. 
This last contribution to the energy must be broken into several portions. It was 
the evaluation of \RR which baffled me, until M. Foincaru’s solution came to 
my help. 

If we imagine the ellipsoid J to be intersected by a family of orthogonal quadrics, 
and if we suppose for the moment that the region R is filled with positive density, 
we may further imagine the matter lying inside any orthogonal tube to be transported 
along the tube, and to be deposited on the surface of J in the form of a concentration 
of positive surface density + 0. The mass of + C is equal to that of + F, but it 
is differently arranged. In the actual system R is filled with negative volume 
density, and we may clearly add to this two equal and opposite surface densities 
+ C and — C on «/. 

Thus the matter lying in the region R may be regarded as consisting of negative 
surface density — C on J, together with a double system, namely negative volume 
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density — R in. conjunction with equal and opposite surface density + C. This 
double system, say D, is therefore C —/A The lost energy \RR may he considered 
as consisting of three parts; first the energy of — C with itself, say ^CC; 
secondly that of D with itself, say ^ DD ; thirdly that of — C with D. This third 
item is obviously equal to — CC + CR , and therefore ^RR is equal to — \CC 
+ CR + iDD. 

It follows that the gravitational lost energy of the pear may be written symboli¬ 
cally in the form 

\JJ - JR + CR- ±CC +' ±I)1X 

In this discussion no attention has as yet been paid to the rotation, but fortunately 
xt happens that the introduction of this consideration actually simplifies the problem, 
for if we suppose \ JJ and JR to mean the lost energies of J with itself and with R 
on the supposition that the mass is rotating with the angular velocity of the critical 
Jacobian, the formulae become much more tractable than would have been the case 
otherwise. 

The inclusion of part of the angular velocity in this portion of the function if, 
only leaves outstanding the excess of the kinetic energy of the pear above the 
kinetic energy, which it would have if it rotated with the angular velocity of the 
critical Jacobian. If co denotes the latter angular velocity, and (co 2 + 8co 2 ) 1 the actual 
angular velocity of the pear; if Aj be the moment of inertia of J , and A r that of R 
considered as filled with positive density, we have 

E = \JJ -JR+ CR - \CC + %DD + i (Aj - A r ) 8 

In this statement I have omitted a term which arises from the displacement of 
the centre of inertia from the centre of the ellipsoid; it is duly considered in the 
paper, but is shown to vanish to the requisite order of approximation (§§ 2, 14). 

The co-ordinates of points are determined by reference to the ellipsoid «/, which 
envelopes the whole pear, and the formula for the internal gravitation of J , inclusive 
of the rotation co, is of a simple character. The size of J is indeterminate, and 
therefore the formulae must involve an arbitrary constant expressive of the size of J ’ 
But the final result E cannot in any way depend on the size of the ellipsoid which 
is chosen as a basis for measurement, and therefore this arbitrary constant must 
ultimately disappear. Hence it is justifiable to treat it as zero from the beginning. 
It appears then that we are justified in using the formula for internal gravity 
throughout the investigation. If the artifice of the enveloping ellipsoid had not 
been adopted, it would have been necessary to take note of the fact that the pear 
is in part protuberant above and in part depressed below the ellipsoid of reference. 
M. Poincare did follow this last plan, and then proceeded to prove the justifiability 
of using the formula for internal gravity throughout, The argument adduced above 
seems, however, sufficient to prove the point. 
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Although the constant expressive of the size of J is put equal to zero—which 
means that the pear is really partly protuberant above the ellipsoid—I have found 
that a considerable amount of mental convenience results from always discussing the 
subject as though the constant were not zero, so that the ellipsoid envelopes the 
pear, and I shall continue to do so here. 

When an ellipsoid is deformed by an harmonic inequality, the volume of the 
deformed body is only equal to that of the ellipsoid to the first order of small 
quantities. In the case of the pear, all the inequalities, excepting the third zonal 
one, are of the second order, and as far as concerns them the volumes of J and 
of the pear are the same. But it is otherwise as regards the third zonal harmonic 
term, and the first task is to find the volume of such an inequality as far as e 1 . 
When this is done we can express the volume of J in terms of that of the pear, 
which is, of course, a constant (§§ 3, 4). 

By aid of ellipsoidal harmonic analysis we may now express the first four terms 
of E in terms of the mass of the pear, and of certain definite integrals which depend 
on the shape of the critical Jacobian ellipsoid (§§ 5, 6, 7). 

The energy J 1)1) presents much more difficulty, and it is especially in this that 
M. Poincare's insight and skill have been shown. The system 1) consists of a layer 
of negative volume density, coated on its outer surface with a layer of surface 
density of equal and opposite mass. 

Two surfaces, infinitely near to one another, coated with equal and opposite surface 
densities, form together a magnetic layer or a layer of doublets. The change of 
potential in crossing such a layer is 4 tt times the magnetic moment at the point of 
crossing, and is independent of the form of surface. To find the difference between 
the potential at two points at a finite distance apart, one being on one side and the 
other on the other side of the layer, we have to add to the preceding difference 
a term equal to the force on either side of the magnetic layer multiplied by the 
distance between the two points. This additional term is small compared with that 
involving the magnetic moment, provided that the distance is small. If the magnetic 
layer coincided with the surface of an ellipsoid the force in question would be exactly 
calculable, and if it lies on the surface of a slightly deformed ellipsoid the force 
remains unchanged by the deformation as a first approximation. 

Thus it follows that it is possible to calculate the difference of potential at two 
points lying on a curve orthogonal to an ellipsoid, when one point is on one side 
and the other on the other side of a magnetic layer residing on a deformation of the 
ellipsoid. Further, if the two points lie on the same side of the magnetic layer the 
term dependent on magnetic lAoment (which would represent the crossing of the 
layer) disappears, and only the term dependent on the force remains. 

Two equal and opposite layers of matter at a finite distance apart may be built up 
from an infinite number of magnetic layers interposed between the two surfaces. 
Hence by the integration of the result for a magnetic layer we may find the change 
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of potential in passing from any one point to any other lying on the same orthogonal 
curve in the neighbourhood of a finite double layer. 

Again, the system D, consisting of — R and + C, may be built up by an infinite 
number of finite double layers. Hence by a second integration we may find the 
difference between the potential of D at any point inside R and the point lying on J 
where the orthogonal curve through the first point cuts the surface of J. 

Finally, it may be proved that the lost energy \DD is equal to half the difference 
of potentials just determined multiplied by the density and integrated throughout 
the region R. The. required expression of this portion of the energy is found to 
consist of two parts, of which one depends on magnetic moment and the other on the 
force (§ 9). The reduction of this part of the energy to calculable forms is not very 
simple ; it is carried out in §§11, 12. 

The calculation of the moment of inertia of the pear is comparatively easy, since it 
only involves those harmonic inequalities of J which are expressible by harmonics of 
the second degree (§ 13). On multiplying the moment of inertia by |-S&b, we obtain 
the last contribution to the expression for E. 

The energy function cannot involve e 3 , since the vanishing of the coefficient of that 
term is the condition whence the critical Jacobian was determined. If /denotes the 
coefficient of any harmonic inequality other than the third zonal one, the part of E 
independent of Sob is found to contain terms in e 3 , c 2 / and ( f )'\ The coefficient of 
Sob consists of a constant term, a term in e 3 and terms in f and /j 3 , where these f’s 
denote the coefficients of the second zonal and sectorial harmonics. This last part 
does not contain the coefficient of any harmonic of odd degree, and in the first part 
the coefficient of the term in e 3 / for all such harmonics is found to vanish. 

The condition for the figure of equilibrium is that the variations of E for variations 
of e 3 and of each / shall vanish. On differentiating E with respect to the / of any 
harmonic of odd degree and equating the result to zero, we see that that / must 
vanish. Hence it follows that the pear cannot involve any odd harmonic excepting 
the third zonal one. Again, the symmetry of the figure negatives the existence of 
any even functions involving sine-functions of the quasi-longitude measured from 
the prime meridian (as I may call it) of symmetry through the axis of rotation, 
ihe same consideration negatives the existence of even functions involving cosine 
functions of odd rank. Accordingly the only functions to be considered are the even 
ones of even rank, involving the cosine functions of the longitude. 

The equation to zero of the variations of E for all the f’s, excepting f, / 3 3 , gives 
at once all those f s in terms of e 3 . The equations to zero of the variations for g 3 , f, 
f* give three equations for the determination of Sob, f, fj as multiples of e 3 . We 
thus have the means of finding the angular velocity and all the f’s in terms of the 
parameter e, which measures the amount of departure of the pear from the critical 
Jacobian ellipsoid (§ 14). 

It seems unnecessary to give here any explanation of the methods adopted for 
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reducing the analytical results to numbers, and it may suffice to say that the task 
proved to be a very laborious one. 

The harmonic terms included in the computation were those of degree 2 and ranks 
0 and 2, of degree 4 and ranks 0, 2, 4, and of degree 6 and ranks 0, 2, 4. The sixth 
sectorial harmonic was omitted because its contribution would certainly prove 
negligible. 

The expression for 8was found in the form of a fraction, of which the denominator 
is determinate and the numerator consists of the sum of an infinite series. Nine 
terms of this series were computed, namely, a constant term and the contribution of 
the eight harmonic terms specified above. I found, in fact, that it would only change 
the numerator by about one-twentieth part of itself, if all the harmonics excepting 
the zonal ones of degrees 2, 4, 6 had been dropped. 

The result shows that the square of the angular velocity of the pear is less than 
that of the critical Jacobian ellipsoid in about the proportion to 1 — -|e 3 to 1. On 
the other hand the angular momentum of the pear is greater than that of the 
ellipsoid in about the proportion of 1 + yy-c 3 to 1. If this last result were based 
on a rigorous summation of the infinite series, it would, in accordance with the 
principle explained above, absolutely prove the stability of the pear. The inclusion 
of the uncomputed residue of the series would undoubtedly tend in the direction 
of reducing the coefficient given in round numbers as yv , and if it were to reduce it 
to a negative quantity, we should conclude that the pear was unstable after all. 
The apparently rapid convergence of the series seemed to render it almost incredible 
that the inclusion of the residue could bring about such a reversal of our conclusion, 
yet I thought it advisable to make a rough estimate of the amount of change which 
would arise from the contribution of* the eighth zonal harmonic. 

The contribution of the sixth zonal harmonic to the series above referred to was 
about ‘00006, and I find that if the contribution of the uncomputed residue should 
amount to ‘00014, the apparent stability of the pear would be just reversed. Now 
my estimate of the contribution of the eighth zonal harmonic to the same sum is 
‘0000008, or only T ygth of the critical amount. 

Since the convergency of the series is obviously very rapid, it is wholly incredible 
that the inclusion of the uncomputed residue could materially alter, much less 
reverse our result. I regard it then as proved, but by something short of an 
absolute algebraic argument, that the pear-shaped figure is stable. 

The numbers obtained in the course of the determination of the stability afford 
the means of giving a second approximation to the form of the pear. The result is 
shown graphically in the figure of § 20, where the largest value of e is adopted 
which seemed to secure a fair degree of approximation in the result. I originally 
called the figure “ pear-shaped,” because M. Poincare’s conjectural sketch in the 

Acta Mathematica ” was very like a pear. In the first approximation, shown in my 
former paper, the resemblance to a pear was not striking, and it needs some imagiua- 
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tion to recognise the pear shape in the second approximation shown here ; but a 
distinctive name is so convenient that we may as well continue to call it by that 
name. 

The effects of the new terms now added are almost entirely concentrated at the 
two ends. All these terms, excepting a very small one arising from the second 
sectorial harmonic, tend to augment the protuberance at the stalk and to fill up the 
depression at the blunt end. It is true that there is a small term, arising from the 
square of the third zonal harmonic, which diminishes the protuberance and increases 
the depression, but this cannot be regarded as a new term, since it only represents 
the effect of the fundamental harmonic carried to the second order of small 
quantities. 

The new zonal harmonics furnish by far the most important contributions. The 
second zonal harmonic denotes that the ellipsoid most nearly resembling the pear is 
longer and less broad than the Jacobian. The largest contribution of all is that 
due to the fourth zonal harmonic, and this may be regarded as the octave of the 
second zonal term. A rough estimate shows that the eighth harmonic, or the double 
octave of the second, is still sensible. The sixth harmonic is the octave of the 
fundamental third zonal harmonic, and is the last of the three important terms. 

The general effect is that the protuberance at the stalk of the pear is much 
increased, and the depression at the other end nearly filled up. Over the greater 
part of the whole surface the depressions and protuberances are less conspicuous 
than they were. The nodal lines where the surface of the pear cuts that of the 
ellipsoid are entirely shifted from their former positions. It did not seem worth 
while to attempt to specify their new positions, because the choice of the ellipsoid to 
which we refer influences the result so largely. The ellipsoid on which these figures 
are constructed is that which is called J in this summary. Its volume is a little less 
than that of the pear, so that the protuberances are a little greater in volume than 
the depressions. 

I think it is hardly too much to say, that in a well-developed “ pear ” the 
Jacobian ellipsoid has nearly regained its primitive figure, but subject to a small 
tidal distortion due to the attraction of a protuberance which it shoots forth at 
one end. I venture to give here a conjectural sketch of a further stage of the 
development. 



If we look at this figure and at those drawn by Mr. Jeans in his striking 
investigation of the parallel changes in the shape of an infinite rotating cylinder 
VOL. CC.—A. 2 S 
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(s'rpra, p. 67), we can hardly fail to be reminded of some such phenomenon as the 
protrusion of a filament of protoplasm from a mass of living matter. 

Notwithstanding the caveat which M. Poincare enters as to the dangers of 
applying these results to heterogeneous masses and to cosmogony, 1 cannot restrain 
myself from joining him in seeing in this almost life-like process a counterpart to at 
least one form of the birth of double stars, planets, and satellites. 


Note. — Erratum in “ Ellipsoidal Harmonic Analysis/’ ‘ Phil. Trans./ A, vol. 197, 
p. 512, line 4 from foot. The first term inside the bracket should be negative. 
The mistake runs on and the same correction should be made in equations 
(62), (63), and (64), and in line 9 on the following page. 



